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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non- SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows; 


_ Multiply _ 

feet 

inches 

kip-foot 

kips  per  square  inch  (ksi) 


_ _ 

0.3048 

0.0254 

1355.818 

6894.757 


To  Obtain 
meters 
meters 

newton-meter 

kilopascals 
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COMPUTER-AIDED.  ELD -VERIFIED 


STRUCTURAL  EVALUATION 

DEVELOPMENT  OF  COMPUTER  MODELING  TECHNIQUES 
FOR  MITER  LOCK  GATES 

PART  I :  INTRODUCTION 

1.  There  are  over  160  commercially  active  locks  (including  over 

200  lock  chambers)  on  the  United  States  inland  and  intracoastal  waterway  sys¬ 
tems.  The  ages  of  these  lock  chambers  range  from  one  year  to  over  150  years. 
Approximately  40  percent  of  the  locks  are  over  50  years  old  and  the  median  age 
of  all  chambers  is  approximately  35  years  (USACE  1988) .  Most  of  the  locks 
include  miter  lock  gates  which  ar  :ither  horizontally  framed  or  vertically 
framed.  Due  to  the  deterioraf  with  age  and  the  large  number  of  miter  gates 
which  are  in  operation,  struct  .i.al  evaluation  is  necessary  to  assess  existing 
strength  and  safety  of  many  of  these  lock  gates.  The  primary  goal  of  the 
project  reported  herein  is  to  develop  an  evaluation  system  incorporating  ana¬ 
lytical  and  experimental  techniques  which  may  be  used  to  assess  the  strength 
and  safety  of  miter  lock  gates. 

2.  Through  the  combined  use  of  modern  computer  software,  digital  data 
acquisition  technology,  and  experimental  testing,  a  simple  realistic  approach 
to  determine  the  structural  behavior  of  miter  lock  gates  is  possible.  Rather 
than  basing  a  structural  evaluation  solely  on  analytical  techniques  or  experi¬ 
mental  measurements,  the  proposed  method  integrates  the  two.  Sophisticated 
finite  element  analytical  models  can  be  developed  to  predict  the  behavior  of 
any  structure.  However,  for  any  analysis  many  assumptions  must  be  used  in  the 
representation  of  structural  geometry  and  material  properties,  and  it  is  not 
possible  to  exactly  represent  the  behavior  of  complex  structures.  Addition¬ 
ally,  sophisticated  experimental  techniques  are  available  to  measure  a  struc¬ 
ture's  response  to  various  loadings,  but  for  a  given  experiment  only  a  few 
selected  points  on  a  structure  can  be  monitored.  By  combining  analytical  and 
experimental  techniques  an  optimum  system  which  incorporates  the  benefits  of 
both  techniques  can  be  developed.  Analytical  and  experimental  results  can  be 
compared,  and  the  analytical  model  can  be  systematically  modified  until  it 
reproduces  the  behavior  observed  under  experimental  conditions. 
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3.  Such  an  integrated  system  has  been  developed  and  successfully  imple¬ 
mented  for  highway  bridges  by  the  University  of  Colorado.  This  work  was 
accomplished  during  the  completion  of  two  research  projects  sponsored  by  the 
Pennsylvania  Department  of  Transportation  (PennDOT)  and  the  Federal  Highway 
Administration  (FHWA)  (Goble,  Schulz,  and  Commander  1990a  and  Goble  et  al . 
1991) .  A  principal  goal  of  these  studies  was  to  develop  a  systematic  approach 
to  evaluating  ambiguous  structural  parameters  so  as  to  improve  the  accuracy  of 
analytical  models.  The  quality  of  the  models  was  evaluated  by  comparing  theo¬ 
retical  structural  responses  with  experimental  (field)  strain  measurements. 

4.  For  the  current  study,  the  modeling  and  analysis  procedures  that 
were  used  for  highway  bridges  have  been  modified  for  use  in  the  evaluation  of 
vertically  and  horizontally  framed  miter  gates.  This  report  describes  the 
methods  used  to  represent  miter  gates  analytically,  the  structural  analysis 
capabilities  of  the  software,  and  the  results  of  two  case  studies  performed  on 
miter  gates  for  which  field  data  had  been  acquired. 
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PART  II;  NUMERICAL  REPRESENTATION  OF  MITER  GATES 


Objectives 

5.  One  of  the  primary  goals  in  the  initial  phase  of  this  project  is  to 
develop  the  analytical  modeling  procedures  required  to  represent  the  struc¬ 
tural  responses  of  miter  gates  subjected  to  normal  service  loads.  The  model¬ 
ing  procedure  should  be  sufficiently  detailed  to  simulate  the  indeterminate 
load  transfer  characteristics  of  miter  gates  yet  be  efficient  and  straightfor¬ 
ward  enough  to  be  used  on  a  routine  basis. 

6.  In  addit_on  to  providing  modeling  and  analysis  procedures,  an 
approach  for  making  modifications  to  the  analytical  model  until  its  load 
response  is  similar  to  that  of  the  actual  structure  is  required.  Quite  often, 
the  numerical  representation  of  an  actual  structure  is  developed  using  a  con¬ 
siderable  amount  of  subjective  judgment.  For  example,  the  structural  geometry 
might  be  simplified  or  a  material  or  structural  property  might  not  be  well 
defined,  particularly  when  representing  a  complex  component  or  a  damaged  ele¬ 
ment.  It  is  desirable  to  quantify  these  obscure  properties  such  that  a  numer¬ 
ical  model  can  provide  load  responses  similar  to  those  of  the  real  structure. 
This  requires  a  method  of  comparing  measured  responses  with  analytical  results 
and  systematically  modifying  the  modeling  parameters  until  the  two  correlate. 
With  this  type  of  iterative  process,  the  final  analytical  model  provides  a 
relatively  accurate  representation  of  the  actual  lock  gate.  The  response  of  a 
structure  under  various  load  conditions  can  then  be  simulated  using  the  ana¬ 
lytical  computer  model,  and  predictions  for  response  can  be  made  for  any  loca¬ 
tion  on  the  structure.  With  this  approach,  the  overall  structural  performance 
can  be  evaluated  based  on  the  correlation  with  experimental  results  obtained 
from  a  few  locations. 


Selection  of  Analysis  Software 

7.  Miter  gates  are  complex  indeterminate  structures,  therefore,  a 
finite  element  approach  is  the  most  logical  choice  for  an  analysis  method. 
Finite  element  procedures  are  relatively  well  understood  and  widely  used  in 
analysis  today. 

8.  Due  to  the  unique  nature  of  the  project  goals,  commercial  software 
packages  that  fulfill  the  requirements  for  both  modeling  and  iterative 
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parameter  modification  do  not  exist.  However,  software  has  been  developed  at 
the  University  of  Colorado  to  address  these  objectives  (Goble,  Schulz,  and 
Commander  1990a  and  Goble  et  al .  1991).  The  computer  program  developed  at  the 
University  of  Colorado,  Structural  Analysis  and  Correlation  (SAC),  meets  most 
of  the  requirements  of  this  project  since  it  has  advanced  modeling  capabili¬ 
ties  and  includes  a  parameter  optimization  algorithm.  To  account  for  geometry 
and  loading  characteristics,  various  features  were  added  to  customize  the 
program  for  miter  gate  modeling  and  analysis. 


Analysis  Capabilities  and  Limitations 


Structural  geometry  and  element  library 

9.  Although  it  is  desirable  to  maintain  as  much  simplicity  as  possible 
in  the  modeling  of  lock  gates,  SAC  includes  three-dimensional  modeling  capa¬ 
bility.  Nodal  displacements  are  represented  with  translations  in  three 
orthogonal  directions  and  rotations  about  three  orthogonal  axes  (six  degrees 
of  freedom  per  node) .  The  overall  structural  geometry  is  defined  within  a 
single  global  Cartesian  coordinate  system,  and  a  local  coordinate  system  is 
included  for  each  individual  element.  Thus,  strains  and  internal  forces  can 
be  represented  with  respect  to  element  orientations  and  there  are  no  limita¬ 
tions  to  the  orientation  of  a  member  within  a  structure. 

10.  In  the  current  version  of  SAC,  structures  can  be  represented  with 
the  following  three  basic  modeling  components: 

a.  Frame  elements.  Girders,  beams,  and  diagonal  members  are 
simulated  with  space  frame  elements.  These  elements  are  one¬ 
dimensional  in  geometry  since  they  are  defined  by  two  nodes; 
however,  they  have  stiffness  properties  to  resist  loads  and 
displacements  in  three  dimensions. 

b.  Plate -membrane  elements.  Skin  plates,  and  for  refined  modeling, 
girder  webs  and  flanges  can  be  modeled  with  the  use  of  rectangu¬ 
lar  membrane -plate  elements.  The  membrane-plate  element  is  a 
hybrid  of  two  separate  quadrilateral  elements:  a  rectangular 
plate  bending  element  and  a  four-node  membrane  element.  The 
plate  bending  element  is  capable  of  modeling  an  out-of-plane 
displacement  and  rotations  about  two  in-plane  axes.  The  mem¬ 
brane  degrees  of  freedom  consist  of  translations  along  the  two 
in-plane  axes.  Because  of  the  small  displacement  assumption, 
the  deflections  are  uncoupled. 

c.  Spring  elements.  Spring  elements  which  may  be  used  to  represent 
elastic  supports  or  semirigid  connections  are  also  available. 
These  elements  have  no  spacial  dimensions  and  are  defined  by 
translational  and/or  rotational  elastic  stiffness  constants. 
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Load  applications 

11.  The  version  of  SAC  developed  for  bridge  studies  was  modified  to 
include  capabilities  to  model  the  primary  loading  characteristics  relevant  to 
miter  gates.  A  strong  effort  has  been  made  to  minimize  input  requirements 
when  implementing  common  load  combinations.  The  loading  applications  cur¬ 
rently  available  include: 

a.  Dead  load.  ''’he  effect  of  gravity  can  easily  be  applied  to  a 
structural  model.  The  material  weight  of  each  element  group  is 
specified  within  a  material  group  input  data  block,  and  the  dead 
load  is  simulated  by  specifying  the  direction  of  gravity.  The 
self  weight  of  each  structural  element  is  automatically  computed 
and  added  to  the  load  vector. 

b.  Water  loading.  Water  pressure  loads  are  simulated  in  the  models 
by  specifying  the  water  elevation  at  the  upstream  and  downstream 
sides  of  the  gate,  the  direction  in  which  depth  is  measured,  and 
the  unit  weight  of  water.  A  constant  pressure  may  also  be 
applied.  Pressure  forces  are  automatically  computed  for  all 
plate  elements  a-iU  are  represented  by  equivalent  nodal  forces 
which  act  normal  to  the  element  plane.  The  magnitude  of  the 
pressure  is  based  on  the  average  depth  of  the  loaded  element. 

c.  Nodal  loads.  Concentrated  forces  and/or  moments  can  be  applied 
at  all  nodal  locations. 

d.  Beam  loads .  Distributed  and/or  concentrated  loads  can  be 
applied  to  specified  frame  elements. 

Loading  capabilities  that  are  likely  candidates  for  future  implementation 
include  temperature  and  prestressing  effects. 

Analysis  results 

12.  Equally  important  to  load  modeling  capabilities  is  the  presentation 
of  analysis  results.  The  number  of  response  types  and  the  clarity  of  the 
presentation  are  vital  to  the  evaluation  process.  SAC  provides  a  complete 
listing  of  analysis  results  in  an  easy  to  read,  page  style  format.  There  are 
print  options  for  each  response  type,  so  only  the  output  desired  will  be  pre¬ 
sented.  The  following  response  types  are  presented  in  an  output  data  file 
with  ASCII  format; 

a.  Input  echo.  To  verify  data  input,  the  geometry,  material  infor¬ 
mation,  boundary  conditions,  and  load  data  are  presented  in  a 
readable  format. 

b.  Nodal  displacements.  As  with  most  finite  element  programs,  the 
displacements  and  rotations  are  tabulated  for  each  node. 

c.  Element  end  forces .  Member  end  forces  are  computed  for  all 
degrees  of  freedom  at  each  element  node  for  both  frame  and  plate 
elements.  End  forces  are  presented  in  relation  to  the  local 
coordinate  system  of  each  individual  element. 
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d.  Strains ,  Strain  can  be  computed  at  any  designated  location  on 
specified  beam  and/or  plate  elements.  These  locations  usually 
correspond  to  the  gage  locations  used  during  field  testing. 
Strains  for  frame  elements  resulting  from  axial  force  P  and 
flexural  force  about  the  major  axis  (y-axis)  of  the  element 
are  determined  by 


£ 


g  ^  ^7 
E  Ely 


(1) 


where  a  is  the  stress,  E  the  modulus  of  elasticity,  y  the  dis¬ 
tance  to  the  neutral  axis,  ly  the  moment  of  inertia  about  the  y 
axis,  and  A  the  cross  sectional  area. 

e.  Support  reactions.  Reaction  forces  are  computed  in  global  coor¬ 
dinates  for  each  node  containing  fixed  degrees  of  freedom. 

f.  Summation  of  applied  loads.  The  summation  of  applied  loads  is 
listed  for  all  six  degrees  of  freedom  as  F^,  Fy,  F^ ,  My,  and 

where  F^  is  the  resultant  force  in  the  global  i  direction 
and  is  the  resultant  moment  about  the  global  i  axis.  These 
results  can  be  used  primarily  as  a  check  for  the  load  input 
data . 

g.  Summation  of  applied  loads  and  reactions.  Modeling  errors  and 
typographical  errors  in  the  input  file  will  lead  to  erroneous 
predictions  such  as  for  structural  instability.  The  summation 
of  the  applied  loads  and  the  reactions  provides  a  check  to 
ensure  that  the  sum  of  forces  is  zero. 

h.  Absolute  and  percent  error.  When  test  data  are  available,  it  is 
often  of  interest  to  obtain  a  comparison  of  results  between  the 
theoretical  and  measured  strains.  The  differences  in  the  strain 
results  are  presented  as  a  summation  of  the  absolute  errors  at 
each  gage  position  (total  strain  difference)  and  as  a  percentage 
error . 

Limitations  of  analysis 

13.  The  main  limitations  to  the  current  version  of  SAC  are  that  it  can 
only  perform  static  analysis  and  that  all  responses  are  assumed  to  be  linear 
and  elastic.  Static  analysis  is  appropriate  for  normal  conditions  since  the 
frequency  of  the  applied  loads,  such  as  hydraulic  loads  due  to  lock  filling  or 
impact  loads,  will  generally  be  much  lower  than  the  natural  frequency  response 
of  miter  lock  gates.  This  is  particularly  true  when  the  gates  are  in  the 
miter  position,  since  this  creates  a  stiff  structural  system  and  thus  has  a 
high  natural  frequency.  The  requirement  of  linear  and  elastic  analysis  is  not 
a  critical  limitation  since  most  load  applications  are  normal  service  loads. 
Data  comparison  and  parameter  evaluation 

14.  The  feature  which  makes  SAC  unique  from  other  structural  analysis 
programs  is  the  systematic  procedure  of  comparing  the  analytical  results  with 
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measured  responses  and  evaluating  parameters  based  on  these  comparisons.  SAC 
includes  an  optimization  procedure  which  is  used  to  minimize  the  absolute 
difference  between  experimentally  obtained  and  theoretical  responses.  Indi¬ 
vidual  structural  parameters  such  as  torsional  stiffness  of  frame  elements  or 
flexural  stiffness  of  damaged  sections  that  are  considered  to  be  unknown  can 
be  optimized  within  user  defined  limits  so  that  the  best  correlation  is 
obtained.  The  optimized  values  can  then  be  incorporated  into  the  model  used 
to  evaluate  the  entire  structure.  The  ability  of  the  optimization  procedure 
to  reduce  the  response  error  is  dependent  primarily  on  the  accuracy  of  the 
geometry  and  assumed  boundary  conditions  of  the  original  model.  If  the  origi¬ 
nal  model  is  basically  correct,  the  error  can  be  minimized  and  the  parameters 
can  be  reasonably  evaluated.  However,  if  the  basic  structural  characteristics 
are  not  well  represented,  a  good  correlation  between  field  and  numerical 
results  will  not  be  obtained.  Hence,  the  optimization  process  acts  as  an 
indicator  of  the  model  quality.  A  more  detailed  description  of  the  optimiza¬ 
tion  procedure  and  algorithm  is  presented  in  Appendix  A. 


Representation  of  Miter  Gates 


15.  In  the  development  of  a  modeling  procedure,  it  was  desired  to  main¬ 
tain  sufficient  simplicity  such  that  the  analysis  program  can  effectively  be 
used  on  a  routine  basis  and  can  be  run  on  a  personal  computer  (PC) .  Although 
there  are  limitations,  modern  PC  hardware  and  software  are  extremely  powerful 
and  can  be  used  to  perform  relatively  detailed  analysis  procedures.  (The  SAC 
program  has  been  developed  to  minimize  memory  requirements  and  maximize  per¬ 
formance  as  described  in  Appendix  B.)  Therefore,  the  primary  concern  is  to 
provide  a  method  which  is  easily  used  on  a  routine  basis.  It  was  desired  to 
develop  a  process  for  which  the  complexity  involved  and  time  required  to  gen¬ 
erate  a  finite  element  model  and  the  associated  input  data  are  minimal.  Addi¬ 
tionally,  the  process  must  be  such  that  the  analysis  will  predict  accurate 
results.  Simplicity  of  a  model  can  be  viewed  in  terms  of  the  structural  geom¬ 
etry  or  in  terms  of  the  theoretical  assumptions  and  calculations  required  to 
compute  cross-section  properties.  There  is  a  direct  compromise  between  geome¬ 
try  and  detail.  The  simpler  the  geometry,  the  more  complex  the  individual 
components  (increased  complexity  in  computation  of  member  properties)  must  be 
if  the  same  degree  of  accuracy  is  to  be  maintained. 
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16.  Based  on  several  case  studies  performed  by  the  authors  and  using 
studies  of  a  previous  CASE  project  (Emkin,  Will,  and  Goodno  1987)  as  an 
example,  several  models  of  various  geometry  configurations  were  selected  for 
vertically  and  horizontally  framed  miter  gates.  The  different  models  vary  in 
complexity,  each  having  unique  advantages  and  limitations  which  are  discussed 
below.  The  model  type  used  in  a  particular  case  is  dependent  on  the  complex¬ 
ity  of  the  actual  structure  and  the  information  sought  from  the  analysis.  All 
of  the  models  outlined  below  can  be  implemented  using  SAC.  One  assumption 
applied  to  all  of  the  models  is  that  symmetry  is  assumed  when  the  gates  are  in 
the  miter  position.  Thus,  only  one  gate  leaf  is  modeled  and  roller  supports 
are  applied  to  the  miter  end  of  the  leaf  to  represent  the  boundary  condition 
imposed  by  the  other  leaf.  The  application  of  roller  supports  is  discussed  in 
further  detail  in  Part  III,  "Model  details." 

Vertically  framed  miter  gate 

17.  Considering  load  transfer,  vertically  framed  miter  gates  are  rela¬ 
tively  simple  compared  to  horizontally  framed  miter  gates.  For  a  mitered 
vertically  framed  gate,  there  are  only  a  few  contact  points  that  resist  water 
and  impact  loads.  There  is  generally  only  one  contact  point  at  the  miter  end 
of  the  leaf  near  the  top  girder,  two  contact  points  at  the  quoin  end,  and  a 
bottom  sill  supporting  the  bottom  girder.  Horizontal  loads  applied  to  the 
vertical  members  are  transferred  directly  to  the  top  girder  and  the  bottom 
sill.  A  vertically  framed  structure  can  be  represented  using  a  grid  model  in 
which  all  of  the  vertical  girders  and  beams  are  simulated  by  frame  elements 
and  are  located  in  a  single  plane  (grid  plane).  In  this  approach  plate- 
membrane  elements  which  are  used  to  simulate  the  skin  plate  or  buckle  plates 
are  located  in  the  plane  with  the  vertical  members.  These  elements  sein/e  to 
resist  horizontal  water  loads  as  well  as  transfer  in-plane  loads.  The  top 
girder  and  diagonal  members  are  offset  from  the  grid  plane  with  rigid  links, 
providing  a  realistic  effect  of  member  eccentricity.  The  horizontal  distances 
(perpendicular  to  the  grid  plane)  between  the  structural  members  and  the  work 
line*  closely  match  the  neutral  axis  location  of  the  members.  Since  the  ver¬ 
tical  members  and  the  buckle  plate  elements  lj*>  in  the  same  plane,  the  stiff¬ 
ening  effects  due  to  composite  action  between  ■  ne  plates  and  the  vertical 


*  The  work  line  is  defined  as  the  line  between  the  contact  points  of  the  top 
girder  which  are  located  at  the  quoin  and  miter  blocks. 
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members  are  accounted  for  in  the  moment  of  inertia  calculations  of  the 
vertical  members. 

18.  The  location  of  the  grid  plane,  which  contains  the  skin  plate  and 
vertical  members,  is  determined  primarily  by  convenience  and  subjective  judg¬ 
ment.  This  plane  should  be  defined  so  as  to  minimize  effort  involved  in 
computing  the  composite  cross-sectional  properties  of  the  vertical  beams  and 
girders.  The  two  recommended  positions  for  the  vertical  plane  would  be  at  the 
skin  plate  or  at  a  location  that  is  assumed  to  be  at  the  average  neutral  axis 
position  of  the  vertical  members. 

19.  Calculation  of  cross-sectional  properties  for  composite  indetermi¬ 
nant  members  involves  subjective  assumptions.  An  initial  assumption  must  be 
made  concerning  the  general  structural  response  and  the  interaction  between 
the  individual  components.  If  it  is  assumed  that  the  vertical  members  will 
act  relatively  independent  of  one  another,  then  the  moment  of  inertia  calcula¬ 
tions  should  be  based  on  the  element  cross  section  including  an  effective 
portion  of  the  skin  plate.  The  assumed  effective  skin  plate  width  can  be 
determined  by  assuming  that  the  effective  portion  is  that  which  is  allowed  for 
unstiffened  elements  under  compression  (USAGE  1984).  When  buckle  plates, 
rather  than  flat  plates,  are  present  their  composite  effects  should  probably 
be  ignored.  Conversely,  if  it  is  assumed  that  the  vertical  members  will  act 
together  as  a  unit,  the  moment  of  inertia  computations  should  account  for  the 
additional  stiffness  imposed  by  the  eccentricity  y  of  the  member  neutral  axis 
relative  to  the  grid  plane.  The  effective  flexural  stiffness  I,  of  a  member 
of  cross-sectional  area  A  can  be  calculated  by 

=  -Tna  ♦  ^7^ 

where  is  the  moment  of  inertia  about  the  member  neutral  axis. 

20.  In  this  case,  the  vertical  plane  should  be  positioned  according  to 
an  average  neutral  axis  calculation.  These  two  methods  of  calculating  cross- 
sectional  properties  should  provide  the  lower  and  upper  bounds  of  an  element's 
flexural  stiffness.  It  is  likely  that  the  true  element  stiffness  will  be 
somewhere  between  the  two  assumed  values  and  can  be  determined  through  corre¬ 
lation  with  field  strain  data  and  optimization. 

21.  One  drawback  to  placing  the  plate  elements  in  the  same  plane  as  the 
vertical  members  is  that  the  membrane  forces  due  to  the  composite  bending 
action  of  the  vertical  members  and  the  skin  plate  are  not  represented  by  the 
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plate  elements.  This  means  that  strains  measured  on  the  skin  plate  cannot  be 
used  for  comparison  with  the  analysis.  In  some  cases,  it  may  be  necessary  to 
produce  a  more  realistic  model  in  which  each  frame  element  is  located  accord¬ 
ing  to  its  calculated  neutral  axis.  Figure  1  illustrates  the  representation 
of  a  vertically  framed  miter  gate  using  a  single  plane  geometry. 

Horizontally  framed 

miter  gate,  simple  grid  model 

22.  A  grid  model  is  a  simple  approach  which  may  also  be  used  to  repre¬ 
sent  the  load  transfer  characteristics  of  a  horizontally  framed  miter  gate. 
Figure  2  illustrates  the  simple  grid  approximation  for  horizontally  framed 
miter  gates.  The  horizontal  girders  and  vertical  diaphragms  are  represented 
by  frame  elements  which  are  located  in  a  single  bent  plane  (grid  plane) .  To 
model  the  eccentricity  between  the  horizontal  girder  neutral  axes  and  the  work 
line,  the  grid  plane  is  offset  near  the  quoin  and  miter  ends  as  shown  in 
Figure  2.  Plate -membrane  elements  representing  the  skin  plate  are  also 
located  in  the  grid  plane.  As  with  the  model  of  a  vertically  framed  gate,  the 
diagonal  member  elements  are  offset  to  there  actual  location  using  rigid 
links,  and  the  most  reasonable  location  of  the  grid  plane  is  at  the  calculated 
average  of  the  horizontal  girder  neutral  axes. 

23,  The  primary  advantage  of  this  approach  is  that  it  is  simple  to 
model  with  respect  to  geometry.  The  disadvantage  is  that  calculations  of 
element  cross-section  properties  become  relatively  difficult,  with  a  large 
degree  of  subjectivity.  Many  assumptions  must  be  made  concerning  the  interac¬ 
tion  between  the  main  structural  elements  and  the  skin  plates  in  computing 
effective  element  stiffnesses.  The  same  disadvantages  inherent  in  the  model 
of  the  vertically  framed  gate  are  present  here.  American  Institute  of  Steel 
Construction  (AISC)  specifications  (AISC  1986,  1989)  provide  limiting  width- 
to- thickness  ratios  for  effective  flange  widths  which  may  be  used  to  approxi¬ 
mate  the  effective  width  provided  by  the  skin  plate.  However,  these  values 
are  design  recommendations  based  on  plate  buckling  requirements  and  are  inher¬ 
ently  conservative.  An  effective  flange  width  equal  to  the  girder  spacing 
provides  a  good  upper  bound  calculation  of  the  girder  stiffness  and  is  likely 
to  be  more  realistic.  As  with  the  vertically  framed  gate,  an  assumption  must 
be  made  concerning  the  interaction  between  the  girders  when  the  neutral  axis 
of  a  girder  cross  section  is  offset  from  the  grid  plane.  If  the  girders  are 
a.ssumed  to  act  independently,  the  moment  of  inertia  about  the  girder  neutral 
axis  should  be  computed.  However,  it  the  girders  are  assumed  to  act  as  a 
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Figure  1.  Representation  of  a  vertically  framed  miter  gate 
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Figure  2.  Simple  grid  model  of  a  horizontally  framed  miter  gate 
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unit,  the  moment  of  inertia  of  girder  cross-section  should  be  calculated  about 
the  reference  plane  according  to  Equation  2. 

24.  Another  consideration  in  determining  cross-section  properties  is 
that  the  cellular  construction  of  the  gate  leaves  tends  to  significantly 
increase  the  torsional  stiffness  of  the  main  structural  members  by  providing 
resistance  against  warping  in  the  girder  and  diaphragm  flanges.  Therefore, 
warping  constants  must  be  implemented  in  the  frame  elements  and  assumptions 
must  be  made  concerning  the  effects  the  skin  plate  has  on  the  warping  torsion. 
If  the  gate  is  represented  in  the  miter  position  or  if  the  majority  of  tor¬ 
sional  strength  of  the  gate  is  assumed  to  be  provided  by  the  diagonal  members, 
then  the  problem  of  torsional  stiffness  in  the  frame  elements  is  of  little 
importance . 

25.  As  with  the  model  of  the  vertically  framed  miter  gate,  calculated 

strains  in  the  skin  plate  are  typically  inaccurate  since  the  membrane  forces 

due  to  composite  bending  are  not  represented  by  this  geometry. 

Horizontally  framed 

miter  gate,  refined  grid  model 

26.  The  refined  grid  model  which  includes  a  more  accurate  representa¬ 
tion  of  structure  geometry  compared  to  the  simple  grid  model  is  illustrated  in 
Figure  3.  It  is  composed  of  the  same  grid  plane  as  that  of  the  simple  grid 
model  with  the  exception  that  the  skin  plate  elements  are  offset  from  the  grid 
with  rigid  links  to  better  represent  the  depth  of  the  gate.  As  with  previ¬ 
ously  described  models,  the  diagonal  member  elements  are  offset  to  their 
actual  location  using  rigid  links.  By  using  this  representation,  the  compos¬ 
ite  behavior  of  the  frame  elements  and  the  skin  plate  is  accounted  for  explic¬ 
itly.  The  effective  increase  in  flexural  member  stiffness  due  to  composite 
action  is  accounted  for  by  the  membrane  action  of  the  plate  elements.  There¬ 
fore,  it  is  not  required  to  include  the  effect  of  the  skin  plate  in  the  frame 
element  moment  of  inertia  calculations  and  the  calculated  strains  in  the  skin 
plate  should  be  more  accurate  than  in  the  simple  grid  approach.  This  approach 
includes  a  more  accurate  representation  of  geometry  than  the  simple  grid 
approach  but  requires  a  significant  increase  in  the  number  of  nodes.  The 
number  of  nodes  in  the  model  is  approximately  doubled  sc  the  amount  of  data 
input  and  the  computer  run  time  are  both  significantly  increased. 

27.  The  calculation  of  girder  and  diaphragm  cross-sectional  properties 
is  simplified  since  the  effects  of  the  skin  plate  do  not  need  to  be  included 
in  the  element  stiffness  computations;  however,  defining  the  location  of  the 
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Figure  3.  Refined  grid  model  for  a  horizontally  framed  miter  gate 
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grid  plane  still  requires  some  assumptions.  Although  several  viable  alterna¬ 
tives  exist,  the  calculated  average  of  girder  neutral  axes  is  the  recommended 
choice  for  the  grid  plane  location.  For  members  with  neutral  axes  eccentric 
to  the  grid  plane  location,  the  eccentricity  should  be  accounted  for  in  the 
flexural  stiffness  computations.  Thus,  the  overall  accuracy  of  the  refined 
grid  model  may  not  be  significantly  improved  when  compared  to  the  simple  grid 
mode  1 . 

28.  As  with  the  simple  grid  model,  warping  torsion  constants  must  be 
implemented  in  the  frame  elements.  Since  the  skin  plates  are  offset  from  the 
frame  elements,  their  effect  on  the  torsional  stiffness  of  the  girders  and 
diaphragms  should  already  be  taken  into  account. 

Horizontally  framed  miter  gate, 
refined  three  dimensional  (3-D)  model 

29.  Another  approach  to  simulating  horizontally  framed  miter  gates  is 
to  model  the  flanges  and  webs  of  each  girder  and  diaphragm  separately.  For 
each  girder  and  diaphragm,  member  flanges  are  represented  by  frame  elements 
and  plate-membratie  elements  are  used  to  simulate  the  web.  The  skin  plate 
elements  are  connected  to  the  flange  nodes  on  the  upstream  face  of  the  gate. 
The  refined  geometry  greatly  simplifies  the  calculation  of  section  properties 
for  the  primary  structural  elements.  The  moment  of  inertia  about  each  axis, 
the  cross-sectional  area,  St.  Venant's  torsion  constant,  and  the  warping  fac¬ 
tor  are  all  implicitly  defined  by  the  geometric  configuration  and  the  elemen¬ 
tary  properties  assigned  to  the  flange  and  web  members.  Figure  4  illustrates 
the  concept  of  modeling  girders  with  frame  and  plate -membrane  elements.  It  is 
apparent  that  the  actual  shape  of  the  structure  is  more  closely  represented 
with  this  approach  than  in  the  two  previously  mentioned  grid  models. 

30.  The  main  disadvantage  of  this  model  is  the  complexity  in  defining 
the  nodal  positions  and  element  connectivity.  The  number  of  nodes  required  is 
approximately  the  same  as  for  the  refined  grid  model,  but  the  number  of  frame 
elements  approximately  doubles  and  the  number  of  plate-membrane  elements 
increases  by  a  factor  of  about  three.  Thus,  user  input,  omputer  memory 
requirements,  and  run  times  increase  significantly.  Implementation  of  this 
modeling  procedure  would  likely  approach  the  memory  limitations  of  most  per- 
.sonal  computers.  It  is  expected  that  memory  requirements  for  most  miter  gate 
models  of  this  type  would  be  between  4  and  8  megabytes. 
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Figure  4.  3-D  representation  of  two  horizontal  girders 


PART  III;  CASE  STUDIES 


31.  To  aid  in  development  of  the  modeling  and  analysis  procedures  as 
applied  to  lock  gates,  two  case  studies  were  performed.  The  first  case  study 
involved  a  vertically  framed  miter  gate.  Experimental  data  were  available  for 
a  barge  impact  test  conducted  on  the  lower  vertically  framed  miter  gate  at  the 
old  Lock  and  Dam  26  north  of  St.  Louis,  MO.  The  test  was  destructive  in  the 
area  of  impact.  However,  strains  measured  at  other  points  on  the  structure 
provided  a  reasonable  basis  for  comparison  since  the  global  behavior  was  elas¬ 
tic  (Chasten  and  Ruf  1991)  . 

32.  A  second  case  study  was  conducted  for  a  horizontally  framed  gate. 

A  single  measurement  of  leaf  dead  load  displacement  was  available  for  a  gate 
leaf  at  the  John  Hollis  Bankhead  Lock  and  Dam  near  Birmingham,  AL  (Emkin, 

Will,  and  Goodno  1987).  This  measurement  was  used  to  indicate  the  suitability 
of  the  modeling  procedures  for  dead  load  response  of  horizontally  framed  miter 
gates.  In  the  following  paragraphs,  each  of  these  case  studies  is  discussed 
in  detail. 


Lock  Gate  26.  Mississippi  River 

33.  The  old  Mississippi  River  Locks  and  Dam  26,  located  just  north  of 
St.  Louis  was  recently  replaced  by  the  Melvin  Price  Lock  and  Dam.  Prior  to 
demolition  of  the  old  locks,  destructive  impact  tests  were  performed  on  the 
downstream  gate  of  the  main  lock  (Goble,  Schulz,  and  Commander  1990b).  The 
gate  which  was  tested  was  a  vertically  framed  miter  gate.  During  the  impact 
tests  the  gate  was  in  the  miter  position  and  no  head  differential  was  present. 
The  impact  load  was  applied  by  a  moving  nine-barge  tow  in  four  separate  occur¬ 
rences.  During  each  test,  strain  data  were  recorded  at  27  locations  on  the 
miter  gate  along  with  the  impact  force,  barge  velocity,  and  barge  accelera¬ 
tion.  Symmetry  of  load  transfer  was  verified  by  nearly  equal  strain  measure¬ 
ments  at  similar  locations  on  both  gate  leaves. 

34.  Data  from  these  tests  provided  a  basis  for  developing  and  evaluat¬ 
ing  analysis  modeling  techniques  for  vertically  framed  miter  gates.  Various 
approaches  in  representing  the  gate  leaf  were  attempted,  starting  with  a  very 
simple  planar  grid  model  in  which  all  elements  were  defined  within  a  single 
vertical  plane.  Modifications  to  the  model  were  implemented  until  a  reason¬ 
able  correlation  between  analytical  and  experimental  results  was  reached.  The 
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process  of  geometry  modification  was  one  of  trial  and  error  and  progressed  as 
the  understanding  of  miter  gate  performance  improved.  In  the  development  of 
the  selected  model,  it  was  apparent  that  the  frame  elements  representing  the 
top  girder  should  be  located  as  close  as  possible  to  the  neutral  axis  of  the 
top  girder.  This  is  necessary  since  the  eccentricity  between  the  line  of 
action  of  the  axial  load  (work  line)  and  the  girder  neutral  axis  results  in  a 
substantial  bending  moment  in  the  top  girder.  The  model  which  was  finally 
selected  consisted  of  the  geometry  discussed  in  Part  II,  "Vertically  framed 
miter  gate." 

35.  Since  the  response  of  the  gate  leaves  in  the  impact  test  was  sym¬ 
metric,  only  one  gate  leaf  was  modeled.  A  concentrated  load  equal  to  one  half 
of  the  measured  impact  load  was  applied  to  the  miter  girder  at  the  point  of 
contact.  Computations  for  the  top  girder  cross-sectional  properties  were 
based  on  the  girder  web  and  flanges  only  (effect  of  skin  plate  was  not 
included) .  The  eccentricity  of  the  top  girder  elements  with  respect  to  the 
work  line  induced  the  appropriate  bending  moment  and  axial  force  responses  in 
the  cop  girder.  A  small  fraction  of  the  top  girder  axial  resultant  force  was 
asstuned  to  be  carried  by  the  skin  plate.  To  accurately  model  this  effect,  it 
was  determined  that  the  grid  plane  (containing  the  buckle  plate  elements  and 
vertical  member  elements)  should  be  located  at  the  center  line  of  th^  buckle 
plates.  However,  this  procedure  may  not  be  the  most  appropriate  for  other 
load  conditions  which  induce  substantial  flexure  in  the  vertical  members. 

Rigid  link  elements  were  used  to  space  the  top  girder  and  the  grid  plane 
appropriately,  and  to  maintain  continuity  between  the  top  girder  and  the  other 
structural  components.  The  miter  and  quoin  contact  points  were  located  at  the 
ends  of  the  top  girder.  Since  the  vertical  beams  and  girders  lie  in  the  same 
plane  as  the  plate  elements,  the  stiffening  effect  of  composite  action  between 
the  plate  and  the  vertical  beams  was  taken  into  account  in  the  moment  of  iner¬ 
tia  calculations  of  the  vertical  members.  However,  the  flexural  stiffness  of 
the  vertical  members  had  little  influence  when  considering  this  loading  situa¬ 
tion  (no  differential  head). 

36.  The  final  model  was  obtained  when  the  correlation  of  the  general 
responses  was  assumed  to  be  accurate,  and  most  results  were  within  20  percent 
of  the  field  measurements.  All  modifications  performed  were  based  strictly  on 
geometry  change! .  No  parameter  optimization  was  performed,  as  the  goal  for 
this  phase  of  the  project  was  to  obtain  general  modeling  characteristics 
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rather  than  structural  parameter  evaluation.  Also,  there  was  no  visible  dete¬ 
rioration  or  any  reason  to  justify  alteration  of  structural  parameters. 

Model  details 

37.  A  description  of  the  analysis  and  modeling  details  is  presented 
here,  and  a  listing  of  the  SAC  input  data  is  provided  in  Appendix  C.  The 
final  model  consisted  of  211  nodal  points  which  translates  into  1266  degrees 
of  freedom  minus  those  eliminated  by  the  boundary  conditions.  There  were 
234  space  frame  elements  and  56  rectangular  plate -membrane  elements.  Twenty 
different  element  groups  were  required  to  define  the  cross-sectional  and  mate¬ 
rial  properties  of  all  the  structural  elements.  All  of  the  material  proper¬ 
ties  for  structural  elements  were  based  on  assuming  36-ksi*  steel.  Rigid  link 
elements  were  used  to  connect  the  top  girder  and  diagonal  elements  to  the  rest 
of  the  structure.  These  elements  were  assigned  extremely  large  axial,  bend¬ 
ing,  and  torsional  stiffnesses  as  compared  to  the  main  structural  elements  and 
are  considered  to  act  in  rigid  manner. 

38.  Since  symmetry  was  ass-..iicd,  only  one  gate  leaf  was  defined  and  the 
action  of  the  other  gate  was  simulated  with  an  effective  inclined  roller  sup¬ 
port  at  the  miter  contact  point.  Inclined  roller  supports  cannot  be  defined 
in  SAC  since  boundary  conditions  are  defined  in  global  coordinates  only. 

Thus,  a  bar  element  aligned  normal  to  the  desired  inclined  roller  plane  was 
used  to  simulate  the  roller  support.  The  bar  element  was  given  axially  rigid 
properties  but  had  no  bending  or  torsional  stiffness.  One  end  of  the  bar 
element  was  fixed  and  the  other  was  located  at  the  miter  contact  point.  To 
assure  that  translation  of  the  contact  point  could  only  occur  along  the 
inclined  plane,  the  bar  element  was  defined  with  a  substantial  length, 
approximately  equal  to  that  of  the  entire  miter  leaf.  To  illustrate  c’-'is 
modeling  concept  further,  a  plane  view  of  the  top  girder  and  the  inclined 
roller  support  is  presented  in  Figure  5.  This  procedure  for  obtaining  symme¬ 
try  can  be  applied  to  horizontally  framed  gates,  regardless  of  the  geometric 
complexity . 

Result  comparisons 

39.  Analysis  results  from  the  final  model  were  compared  with  field  test 
results  at  several  locations  on  the  lock  gate.  Experimental  strain 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  4. 
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Center  Line 
of  Lock 


Figure  5.  Representation  of  roller  supports  to  model  symmetry 


measurements  were  used  to  calculate  the  axial  force  and  flexural  moment  for 
specific  element  cross  sections.  The  emphasis  on  the  comparisons  was  based  on 
responses  of  the  top  girder  since  it  was  responsible  for  transferring  the 
majority  of  the  impact  load  to  the  lock  walls.  The  experimental  results  and 
corresponding  analytical  results  are  presented  in  Tables  1-4.  The  results  are 
presented  for  two  values  of  impact  magnitude  for  each  of  the  four  events. 

Table  1  includes  strain  responses,  corresponding  to  locations  shown  by  Fig¬ 
ure  6,  for  the  middle  of  the  top  girder.  The  strain  comparisons  on  the  web  at 
this  section  may  not  be  very  accurate  as  there  was  a  large  hole  in  the  web  of 
the  girder  near  the  instrumentation.  A  comparison  of  axial  force  and  bending 
moment  at  the  center  of  the  top  girder  is  provided  in  Table  2.  Tables  3  and  4 
contain  the  axial  force  and  moment  comparisons  near  the  miter  and  quoin  ends 
of  the  top  girder.  Strain  measurement  devices  (transducers)  were  placed  at 
approximately  the  quarter  points  of  the  top  girder  as  illustrated  in  Figure  6. 


QUOIN 

GIRDER 


TO 


TD 


MITER 

GIRDER 


BOTTOM 


r 

BOTTOM 


KEY: 

Tronsolwcer 


Figure  6.  Transducer  locations  on  top  girder 


Case  study  conclusions 

40.  The  primary  conclusion  that  was  drawn  from  this  study  is  the  basic 
analysis  and  modeling  procedures  are  appropriate  for  obtaining  reasonably 
accurate  predictions.  Based  on  the  tabulated  results,  it  is  apparent  that 
most  of  the  calculated  responses  are  within  20  percent  of  the  measured 
results.  Considering  the  assumptions  involved  in  modeling  and  in  calculating 
experimental  force  resultants,  this  can  be  considered  a  good  correlation.  The 
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Table  1 

Strains  at  Tod  Girder  (Center  Cross  Section) 


Strain 

( in. /in. 

X  10'®') 

■■1 

mu 

Impact 

Impact 

Load 

H 

K 

L 

M 

G 

Number 

nSISI 

F* 

A* 

F 

A 

F 

A 

F 

A 

F 

A 

1  (max) 

442.9 

-386 

-302 

-109 

-201 

-99 

-109 

-76 

-53 

-8 

-1 

(=75%) 

350.0 

-313 

-239 

-81 

-159 

-67 

-86 

-48 

-42 

-3 

-0.8 

2  (max) 

443.6 

-386 

-303 

-122 

-201 

-105 

-109 

-65 

-53 

15 

-1 

(=75%) 

344.2 

-293 

-235 

-83 

-156 

-70 

-85 

-43 

-41 

5 

-0.8 

3  (max) 

604.7 

-587 

-412 

-144 

-274 

-141 

-149 

-96 

-72 

-16 

-1 

(=75%) 

453.2 

-347 

-309 

-92 

-206 

-82 

-112 

-61 

-54 

-5 

-1 

*  F  =  field,  A  =  analysis. 


Table  2 

Axial  Force  and  Flexural  Moment.  Tod  Girder  Center 


Impact 


Impact 

Load 

Axial  Force 

(kins) 

Bendin£  Moment 

(kin- ft) 

Number 

Field 

Analysis 

%  Error 

Field 

Analysis 

%  Error 

1  (max) 

442.9 

-432.2 

-371.52 

-14.04 

1173.0 

1019.9 

-13.05 

(=75%) 

350.0 

-331.6 

-293.6 

-11.46 

968.3 

806.0 

-16.76 

2  (max) 

443.6 

-396.7 

-372.1 

-6.20 

1266.3 

1021.5 

-19.33 

(=75%) 

344.2 

-302.2 

-288.7 

-4.47 

937.9 

792.6 

-15.49 

3  (max) 

604.7 

-640.9 

-507.2 

-20.86 

1771.4 

1392.5 

-21.39 

tt 

LTI 

453.2 

-377.1 

-380.2 

0.82 

1061.6 

1043.6 

-1.70 
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Table  3 


Axial  Force  and  Flexural  Moment,  Top  Girder  Miter  End 


Impact 


Impact 

Number 

Load 

(kios) 

Axial  Force  (kios) 

Bendine  Moment 

Ckin- ft) 

Field 

Analysis 

%  Error 

Field 

Analysis 

%  Error 

1  (max) 

442.9 

-509.0 

-402.9 

-20.84 

1160.0 

1114.9 

-3.89 

(=75%) 

350.0 

-377.5 

-318.4 

-15.66 

914.7 

881.1 

-3.67 

2  (max) 

443.6 

-379.8 

-403.5 

-6.24 

1382.0 

1116.7 

-19.20 

(=75%) 

344.2 

-360.7 

-313.1 

-13.20 

1145.5 

866.4 

-24.36 

3  (max) 

604.7 

-1621.4 

-550.1 

-66.07 

3989.9 

1522.2 

-61.85 

(=75%) 

453.2 

-1399.9 

-412.3 

70.55 

3580.7 

1140.8 

-68.14 

Table  k 

Axial  Force  and  Flexural  Moment.  Top  Girder  Quoin  End 


Impact 


Impact 

Load 

Axial 

Force  (kips) 

Bendine  Moment 

(kio-ft) 

Niimber 

(kips) 

Field 

Analysis 

%  Error 

Field 

Analysis 

%  Error 

1  (max) 

442.9 

-446.4 

-426.6 

-4.44 

1123.0 

919.2 

-18.15 

(=75%) 

350.0 

-363.0 

-337.1 

-7.13 

873.2 

726.4 

-16.81 

2  (max) 

443.6 

-414.2 

-427.3 

3.16 

1177.5 

920.7 

-21.81 

(=75%) 

344.2 

-332.7 

-331.5 

-0.36 

829.7 

714.4 

-13.90 

3  (max) 

604.7 

-705.1 

-582.5 

-17.39 

1743.9 

1255.0 

-28.03 

(=75%) 

453.2 

-428.3 

-436.8 

1.98 

1039.7 

<540.6 

-9.53 
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results  were  better  at  locations  away  from  the  miter  end  since  there  was 
yielding  near  the  impact  point  (see  Table  3,  impact  3). 

41.  It  was  previously  noted  that  the  final  model  was  obtained  by  modi¬ 
fying  the  model  geometry  to  obtain  the  desired  accuracy.  All  section  parame¬ 
ters  were  computed  using  basic  principles  of  mechanics  of  materials  and  no 
parameter  optimization  was  performed.  It  is  assumed  that  these  parameters  are 
reasonably  accurate  since  the  lock  gate  showed  no  signs  of  significant 
deterioration.  One  structural  parameter,  the  torsional  resistance  of  the 
frame  elements,  may  not  be  highly  accurate  and  may  require  optimization  when 
other  loading  conditions  such  as  dead  load,  water  loading,  or  gate  operation 
are  applied.  Beam  torsion  constants  were  obtained  from  the  basic  equation  of 
nonrestrained  thin-walled  open  sections  (S  hitj^/3) .  For  loading  situations 
when  element  torsion  is  significant,  it  is  likely  that  restrained  warping 
torsion  will  need  to  be  accounted  for.  Lock  gate  geometry  is  such  that  girder 
and  diaphragm  flanges  are  usually  restrained  against  warping  which 
significantly  increases  torsional  stiffness  of  the  gate.  Restrained  warping 
torsion  can  be  modeled  by  either  applying  warping  constant  terms  in  the  tor¬ 
sional  stiffness  calculations  of  the  frame  elements,  or  by  using  a  refined  3-D 
geometry  in  which  the  flange  warping  is  simulated  by  the  additional  degrees  of 
freedom . 


John  Hollis  Bankhead  Lock  and  Dam 


42.  The  John  Hollis  Bankhead  Lower  Miter  Gate  on  the  Black  Warrior 
River  in  Alabama  was  the  subject  of  a  finite  element  study  on  horizontally 
framed  miter  gates  (Emkin,  Will,  and  Goodno  1987).  Reports  generated  from 
this  study  were  reviewed  to  obtain  general  modeling  procedures  and  were  used 
to  provide  a  basis  for  analysis  evaluation.  A  single  field  measurement  was 
obtained  from  this  study  and  was  used  as  the  validation  criteria  for  the 
adopted  modeling  techniques.  An  out-of-plumb  displacement  of  the  gate  leaf 
due  to  gravity  loading  only  was  taken  prior  to  the  attachment  of  the  diagonal 
braces.  It  was  assumed  that  it  was  not  a  highly  accurate  measurement  as 
residual  and  fabrication  stresses  alone  could  account  for  a  considerable 
amount  of  warping  in  the  gate  leaf.  It  was  desirable,  however,  to  obtain 
results  in  the  same  order  of  magnitude. 

43.  Since  the  dead  load  acts  at  the  structure  center  of  gravity  which 
is  eccentric  from  the  structure  shear  center  (indicated  by  the  out-of-plumb 
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displacement) ,  the  refined  grid  model  was  selected  as  the  appropriate  geometry 
configuration.  The  simple  grid  model  was  not  capable  of  accurately  represent¬ 
ing  the  eccentricity  of  the  dead  load,  and  the  3-D  model  was  considered  too 
involved  to  utilize  without  the  aid  of  automatic  mesh  generation  and  was 
therefore  not  attempted. 

Model  development 

44.  Miter  gate  details  were  obtained  from  contract  drawings  and  a 
numerical  model  of  one  leaf  was  developed  using  the  refined  grid  approach 
outlined  in  Part  II  (Figure  3) .  The  girders  and  diaphragms  were  simulated  by 
frame  elements  and  the  skin  plate  was  simulated  using  rectangular  plate- 
membrane  elements.  To  account  for  the  eccentricity  between  the  skin  plate  and 
neutral  axis  of  the  girders  and  diaphragms,  rigid  links  were  utilized  to  off¬ 
set  the  skin  plate  in  a  realistic  manner.  The  girder  and  diaphragm  frame 
elements  were  located  at  the  position  of  the  center  of  area  of  the  girders. 
This  location  was  selected  so  the  load  due  to  structure  weight  would  be  posi¬ 
tioned  appropriately.  Boundary  conditions  were  applied  at  the  pintle  and 
gudgeon  pin  locations  and  an  additional  support  was  added  to  the  bottom  girder 
simply  to  maintain  a  stable  structure  (i.e.  keep  the  gate  leaf  from  swinging). 
Gravity  loads  were  automatically  calculated  and  applied  by  providing  the  unit 
weight  of  steel. 

45.  The  out-of-plumb  field  measurement  (Emkin,  Will,  and  Goodno  1987) 
was  the  horizontal  displacement  of  the  top  miter  corner  relative  to  the  bottom 
miter  corner  and  was  equal  to  4.75  in.  Results  from  the  first  analysis  run 
indicated  an  out-of-plumb  displacement  of  44.7  in.  which  is  in  error  by  a 
magnitvide  of  about  10. 

46.  It  was  determined  that  much  of  the  error  was  due  to  the  extremely 
low  torsional  stiffness  values  given  to  the  frame  elements.  The  torsion  con¬ 
stants  were  computed  for  thin-walled  open  sections  subject  to  St.  Venant's 
uniform  torsion,  in  which  all  of  the  torsion  is  resisted  by  the  shear  stress 
in  the  cross-section  walls  and  resistance  due  to  warping  is  neglected.  How¬ 
ever,  the  cellular  construction  of  the  gate  leaf  induces  a  considerable  amount 
of  resistance  to  flange  warping.  The  consideration  of  restrained  warping  in 
the  flanges  increases  the  overall  torsional  stiffness  of  the  horizontal  gird¬ 
ers  by  a  factor  of  about  20.  Therefore,  the  model  was  modified  to  incorporate 
warping  constants  of  element  cross  sections  in  the  torsional  stiffness  terms. 

A  discussion  of  warping  torsion  as  it  pertains  to  the  analysis  is  presented  in 
Appendix  D. 
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47.  The  warping  constants  were  computed  for  each  frame  element  cross 
section  and  entered  in  the  input  file.  The  new  analysis  resulted  in  an  out- 
of -plumb  displacement  prediction  of  18.6  in.  which  was  a  considerable  improve¬ 
ment  over  the  initial  results.  Implementing  the  restrained  warping  torsion 
terms  in  the  model  significantly  increased  the  torsional  stiffness  of  the  gate 
leaf  and  reduced  the  error  to  a  magnitude  of  4,  but  this  was  still  regarded  as 
unacceptable . 

48.  Further  examination  revealed  a  modeling  deficiency  in  the  rectangu¬ 
lar  plate-membrane  elements.  As  discussed  in  Part  II,  "Structural  geometry 
and  element  library,"  the  hybrid  plate-membrane  elements  are  composed  of  two 
separate  elements;  (a)  a  rectangular  Kirchhoff  plate  element  that  resists  an 
out -of -plane  displacement  and  rotations  about  the  two  in-plane  axes,  and  (b)  a 
plane  strain  quadrilateral  membrane  element  that  resists  displacements  in  the 
two  in-plane  axes.  The  resistance  of  each  of  the  components  is  uncoupled  due 
to  the  small  displacement  assumption.  The  combined  stiffness  terms  in  the 
plate -membrane  element  account  for  only  five  degrees  of  freedom  at  a  given 
node  as  shown  by  Figure  7(A),  and  the  plate -membrane  element  does  not  include 
a  stiffness  term  for  the  sixth  degree  of  freedom  (rotation  about  the  out-uf- 
plane  z-axis).  This  degree  of  freedom  corresponds  to  that  of  the  bending 
stiffness  about  the  weak  axis  of  the  attached  frame  elements  which  influences 
the  overall  gate  leaf  torsional  resistance.  Hence,  by  neglecting  the  sixth 
degree  of  freedom  in  the  plate -membrane  element,  the  numerical  representation 
of  the  gate  leaf  is  too  flexible  when  torsional  resistance  is  a  factor.  Fig¬ 
ure  7(A)  illustrates  the  degrees  of  freedom  resisted  by  the  plate-membrane 
element,  and  the  modeling  deficiency  as  it  applies  to  this  application  is 
shown  in  Figure  7(B). 

49.  To  remedy  this  situation,  two  possible  solutions  exist.  The  first 
is  the  development  and  implementation  of  a  new  plate  element  which  contains 
the  stiffness  terms  for  the  sixth  degree  of  freedom.  The  second  approach 
would  be  to  refine  the  mesh  (reduce  size  of  plate-membrane  elements)  as  shown 
in  Figure  8  so  that  the  in-plane  resistance  of  the  plate  elements  would  be 
effective  in  resisting  weak  axis  bending  of  the  attached  frame  elements. 

These  options  were  not  considered  to  be  feasible  within  the  scope  of  this 
project  and  since  the  deficiency  was  not  deemed  relevant  to  the  majority  of 
load  conditions  they  were  not  pursued. 

50.  Adding  the  sixth  degree  of  freedom  stiffness  terms  to  the  skin 
plate  elements  would  effectively  increase  the  lateral  (weak  axis)  stiffness  of 


30 


Figure  7.  (A)  Degrees  of  freedom  associated  with  plate -membrane  element 

(B)  Modeling  deficiency  due  to  missing  rotational  stiffness  terms 


—  fnccJ  f r  one  elenent 


Figure  8.  Mesh  refinement  approach  to  modeling  in-plane  bending 

the  elements  representing  the  horizontal  girders,  diaphragms,  and  the  quoin 
and  miter  girders.  Therefore,  it  was  of  interest  to  determine  the  effect  of 
increasing  the  lateral  stiffness  of  the  frame  elements.  Although  the  proce¬ 
dure  may  be  purely  academic,  it  was  assumed  that  an  effective  flange  width  of 
U  ft  (spacing  between  the  lower  horizontal  girders)  could  be  contributed  by 
the  skin  plate.  The  moments  of  inertia  about  the  weak  axes  were  then 
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recalculated  using  this  assumption.  Reanalysis  using  the  modified  input  data 
resulted  in  an  out-of-plumb  displacement  of  3.1  in.,  which  was  well  within  the 
desirable  error  limit. 

51.  It  may  be  argued  that  some  of  the  skin  plate  effects  were  now  being 
applied  to  the  gate  leaf  twice,  increasing  the  gate  stiffness  by  an  unrealis¬ 
tic  amount.  As  a  means  of  determining  the  effect  that  the  skin  plate  elements 
have  on  the  torsional  stiffness  of  the  gate  leaf,  the  modulus  of  elasticity  in 
the  skin  plate  elements  was  set  to  zero  and  the  analysis  was  again  performed. 
The  result  was  the  plate -membrane  elements  had  very  little  effect,  as  the  new 
displacement  result  was  evaluated  at  4.9  in.  Compared  with  the  field  measure¬ 
ment  of  4.7  in.,  this  was  considered  to  be  adequate.  Even  though  the  lateral 
stiffness  properties  applied  to  the  model  cannot  be  considered  accurate,  the 
general  effect  of  the  modeling  deficiency  in  the  plate -membrane  elements  was 
verified. 

Case  study  conclusions 

52.  Evaluation  using  the  out-of -plumb  displacement  (measured  with  diag¬ 
onals  not  attached)  is  generally  not  relevant  to  miter  gate  performance  as  the 
majority  of  torsional  stiffness  is  provided  by  the  diagonals  and  the  torsional 
rigidity  is  a  minor  consideration  when  the  gate  is  in  the  miter  position. 

Dead  load  deflections  will  generally  not  be  obtainable  for  existing  structures 
because  the  diagonals  are  already  in  place,  so  the  use  of  the  refined  grid 
model  may  not  be  necessary.  Since  most  loading  is  applied  normal  to  the  plane 
of  the  gate,  this  measurement  should  not  be  the  sole  means  of  evaluating  the 
modeling  procedures,  yet  it  does  provide  some  insight  as  to  the  suitability  of 
the  numerical  model.  An  extremely  good  correlation  was  achieved  with  the 
refined  grid  approach,  although  some  of  the  obtained  cross-sectional  proper¬ 
ties  cannot  be  completely  rationalized.  If  it  is  determined  that  this  type  of 
load  response  is  important,  then  the  3-D  model  would  be  the  best  approach 
since  the  independent  modeling  of  the  beam  flanges  more  realistically  depicts 
the  warping  torsion  problem.  A  new  membrane  element,  capable  of  resisting 
rotation  about  the  out-of-plane  axis  (sixth  degree  of  freedom),  would  also 
need  to  be  developed  and  implemented  to  obtain  a  realistic  overall  torsional 
stiffness . 
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PART  IV:  GENERAL  CONCLUSIONS 


53.  The  approach  of  modeling  and  modifying  the  structural  parameters 
has  been  tested  with  a  high  level  of  success  on  highway  bridges  (Goble, 

Schulz,  and  Commander  1990a  and  Goble  et  al.  1991).  Based  on  these  initial 
case  studies,  it  is  apparent  that  it  is  applicable  to  miter  gates  as  well. 

The  level  of  field  verification  to  date  is  limited  and  a  considerable  amount 
of  experimental  testing  is  still  required.  The  tools  developed  so  far  are 
sufficiently  refined  to  proceed  with  prospective  case  studies;  however,  it  is 
assumed  that  modifications  to  modeling  procedures  may  be  necessary,  pending 
comparisons  with  future  field  tests.  The  analysis  program  can  be  modified  as 
required  to  implement  new  modeling  capabilities. 

54.  Since  field  verification  is  very  limited,  a  clear  decision  regard¬ 
ing  the  best  model! procedures  cannot  be  made  at  this  time.  It  is  proposed 
that  specific  mo  .e'  geometries  will  be  selected  based  on  correlations  made 
with  field  te? rs  which  are  scheduled  for  future  phases  of  this  project.  Con¬ 
clusions  made  at  this  time  are  mainly  qualitative  and  subject  to  change. 

Vertically  Framed  Miter  Gates 

55.  Based  on  this  initial  study,  the  only  modeling  procedure  presented 
for  vertically  framed  miter  gates  is  that  of  the  simple  model.  However,  as  in 
the  case  for  horizontally  framed  miter  gates,  a  3-D  model  may  be  required  for 
detailed  analyses  for  certain  loading  conditions  such  as  those  which  cause 
gate  leaf  twist. 


Horizontally  Framed  Miter  Gates 

56.  The  refined  grid  model  was  selected  for  the  horizontally  framed 
gate  leaf  modeled  in  the  second  case  study.  This  was  due  to  the  nature  of  the 
gravity  load  response  and  the  single  available  displacement  measurement. 
Typically,  this  load  response  will  be  of  little  consequence  so  that  the  use  of 
the  refined  grid  model  may  not  be  the  best  alternative.  When  considering 
model  generation,  the  refined  grid  model  is  considerably  more  complicated  than 
the  simple  grid  model  and  there  is  little  overall  benefit.  For  most  applica¬ 
tions,  the  major  advantages  of  the  refined  grid  are  that  the  theoretical 
strains  in  the  plate  will  be  more  realistic  and  the  composite  action  of  the 
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skin  plate  does  not  need  to  be  considered  in  the  calculation  of  the  horizontal 
girder  properties.  The  same  assumptions  regarding  the  location  of  the  frame 
elements  must  be  applied  in  both  the  simple  and  refined  grid  models,  so  a 
significant  improvement  in  accuracy  may  not  be  achieved  by  using  the  refined 
grid  model.  Requiring  a  fraction  of  the  modeling  effort  and  computer  run 
time,  the  simple  grid  geometry  is  generally  the  preferred  modeling  approach. 
The  primary  disadvantage  of  the  simple  approach  is  in  the  computation  of 
effective  element  stiffness  values  as  discussed  in  Part  II,  "Vertically  framed 
miter  gate."  When  a  more  sophisticated  representation  is  required  such  as 
prediction  of  gate  twist,  it  would  then  be  beneficial  to  implement  the  3-D 
geometry  since  it  is  significantly  more  realistic  than  either  of  the  grid 
models . 
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PART  V;  FUTURE  IMPLEMENTATIONS 


57.  Regardless  of  the  model  geometry  selected,  an  interactive  prepro¬ 
cessor  which  aids  in  rapid  development  and  modification  of  miter  gate  finite 
element  meshes  is  required  to  make  routine  analysis  feasible.  Graphical 
representation  of  structure  geometry,  load  conditions,  and  structural 
responses  would  reduce  the  number  of  user  input  errors  and  aid  in  the  detec¬ 
tion  of  gross  errors  in  geometry  and  boundary  conditions.  A  visual  display  of 
measured  and  calculated  responses  would  also  provide  intuitive  insight  to  lock 
gate  responses.  It  is  intended  to  develop  software  with  these  capabilities 
during  future  phases  of  this  project. 

58.  The  comparisons  of  experimental  and  analytical  data  have  thus  far 
indicated  that  SAC  can  provide  accurate  representation  of  miter  gates.  How¬ 
ever,  the  field  data  at  this  stage  are  rather  limited.  To  fully  verify  the 
modeling  and  analysis  procedures,  data  comparisons  for  a  variety  of  lock  gates 
and  load  configurations  are  recommended. 

59.  Through  the  John  Hollis  Bankhead  Lock  Gate  study,  it  was  observed 
that  a  modeling  deficiency  exists  in  the  plate -membrane  elements  currently 
utilized  by  the  analysis  software.  This  deficiency  is  assumed  to  be  a  minor 
consideration,  generally;  however  upon  further  field  verification,  the  imple¬ 
mentation  of  a  refined  membrane  element  may  be  substantiated. 

60.  Due  to  various  eccentricities  and  member  neutral  axis  locations, 
determining  flexural  stiffness  constants  for  the  horizontal  or  vertical  gird¬ 
ers  and  selecting  appropriate  locations  for  the  grid  plane  is  a  common  problem 
in  the  development  of  each  grid  model.  Using  the  current  models,  elements  for 
members  which  actually  have  different  locations  must  be  located  in  a  common 
plane.  This  is  required  to  make  generation  of  the  models  feasible  and  to 
maintain  continuity  at  nodal  points  with  multiple  elements  attached.  A 
relatively  simple  approach  to  eliminating  this  problem  would  be  to  modify  the 
analysis  program  so  that  eccentricity  terms  are  included  in  the  individual 
element  stiffness  matrices.  The  resulting  benefit  is  that  a  common  reference 
plane  could  be  used  to  define  the  entire  structure,  yet  each  element  could 
have  a  unique  location.  This  would  essentially  provide  the  same  effect  as 
using  rigid  links  to  connect  an  element  to  the  plane,  but  additional  nodes  and 
elements  would  not  be  required.  Implementation  of  the  eccentricity  terms 
would  be  relatively  simple  since  they  can  be  applied  through  transformation 
matrices.  An  example  of  this  type  of  element  is  shown  in  Figure  9. 
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Figure  9.  Frame  element  with  built-in  eccentricity  term 


36 


REFERENCES 


American  Institute  of  Steel  Construction.  1986.  Manual  of  Steel  Construc¬ 
tion,  Load  and  Resistance  Factor  Design  (LRFD) ,  1st  ed. ,  Chicago,  IL. 

_ .  1989.  Manual  of  Steel  Construction ,  Allowable  Stress  Design 

(ASD) .  9th  ed. ,  Chicago,  IL. 

Chasten,  C.  P.,  and  Ruf,  T.  1991.  "Miter  Gate  Barge  Impact  Testing,  Locks 
and  Dam  26,  Mississippi  River,"  Proceedings  of  the  1991  Corps  of  Engineers 
Structural  Engineering  Conference .  (Proceedings  in  preparation) 

Emkin,  L.  Z. ,  Will,  K.  M. ,  and  Goodno,  B.  J.  1987.  Finite  Element  Studies  of 
a  Horizontally  Framed  Miter  Gate,  Technical  Report  ITL-87-4,  US  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  MS. 

Frangopol,  D.M.,  and  Klisinski,  M.  1989.  "Material  Behavior  and  Optimum 
Design  of  Structural  Systems,"  Journal  of  Structural  Engineering,  ASCE , 

Vol  115,  No.  5,  p  1054. 

Goble,  G.  G.,  Schulz,  J.  X.,  Commander,  B.  C.  1990a.  Simple  Load  Capacity 
Tests  for  Bridges  to  Determine  Safe  Posting  Levels:  Final  Report,  submitted  to 
the  Pennsylvania  Department  of  Transportation  by  the  Department  of  Civil  and 
Architectural  Engineering,  University  of  Colorado,  Boulder,  CO. 

_ .  1990b.  Lock  and  Dam  tt26  Field  Test  Report  for  the  Army  Corps  of 

Engineers,  submitted  by  Bridge  Diagnostics,  Inc.,  Boulder,  CO. 

Goble,  G.  G.  ,  Schulz,  J.  X.,  Commander,  B.  C. ,  Burgess,  C.  U. ,  Dow,  J.  0.,  and 
Frangopol,  D.  M.  1991.  Load  Prediction  and  Structural  Response:  Final 
Report,  submitted  to  the  Federal  Highway  Administration  by  the  Department  of 
Civil  and  Architectural  Engineering,  University  of  Colorado,  Boulder,  CO. 

Robson,  B.  1990.  "Predicting  Structural  Response  in  Highway  Bridges,"  M.S. 
Thesis,  University  of  Colorado,  Boulder,  CO. 

Rosenbrock,  H.  H.  1960.  "An  Automated  Method  for  Finding  the  Greatest  or 
Least  Value  of  a  Function,"  Computer  Journal,  Vol  3,  p  175. 

US  Army  Corps  of  Engineers.  1984.  Engineering  and  Design,  Lock  Gates  and 
Operating  Equipment ,  Engineer  Manual  1110-2-2703,  Washington  D.C. 

_ .  1988.  The  1988  Inland  Waterway  Review,  Institute  for  Water 

Resources,  Ft.  Belvoir,  VA. 


37 


APPENDIX  A:  PARAMETER  OPTIMIZATION 


1.  The  unique  feature  of  SAC  is  its  ability  to  compare  certain  analysis 
results  with  a  corresponding  set  of  data  and  to  modify  structural  parameters 
so  as  to  obtain  the  best  correlation.  This  feature  was  developed  to  provide  a 
means  of  evaluating  various  parameters  that  cannot  be  well  defined.  Parameter 
modification  is  performed  using  a  constrained  optimization  technique,  in  which 
the  user  can  specify  what  structural  parameters  are  to  be  optimized.  The 
process  is  considered  constrained  since  the  optimized  values  must  fall  within 
user  specified  limits.  Parameters  which  can  be  optimized  are  the  material 
stiffnesses  and  cross-sectional  properties  for  each  user  defined  element 
group . 

2 .  Variable  optimization  is  performed  by  an  iterative  process  in  which 
the  goal  is  to  minimize  a  single  objective  function.  The  objective  function 
is  defined  as  the  scun  of  the  absolute  differences  between  the  analysis  results 
and  the  corresponding  test  data.  The  structural  response  used  to  calculate 
the  objective  function  is  strain  at  selected  experimental  gage  locations.  The 
absolute  value  of  the  strain  difference  is  computed  at  each  gage  location  and 
then  summed  to  form  the  objective  function.  Thus  any  differences  in  the  theo¬ 
retical  and  experimental  strain  values  will  always  result  as  an  increase  in 
the  objective  function.  The  goal  of  the  optimization  is  to  then  minimize  this 
error  function.  Once  the  analysis  is  run  and  the  error  value  is  defined,  the 
specified  parameters  are  systematically  altered  and  the  analysis  is  run  again. 
The  new  error  is  compared  with  the  previous  value  and  it  is  determined  if  the 
correlation  improved.  This  iterative  process  is  continued  until  an  acceptable 
correlation  (i.e.  an  error  of  less  than  one  percent  of  the  total  absolute 
measured  strain)  is  obtained  or  until  a  maximum  number  of  iterations  are 
reached.  If  an  acceptable  correlation  cannot  be  achieved,  this  indicates  that 
the  model  is  not  an  appropriate  representation  of  the  actual  structure. 

3.  The  optimization  process  used  in  SAC  is  a  "constrained,  direct 
search  method  for  a  single-objective  minimization  algorithm."  This  algorithm 
was  developed  by  Frangopol  and  Klisinski  (1989)*  and  is  an  extension  of  the 
method  of  unconstrained  minimization  by  Rosenbrock  (1960) .  The  integration  of 


*  See  REFERENCES  at  end  of  main  text. 
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this  process  with  an  analysis  program  was  first  performed  during  a  study  spon¬ 
sored  by  the  Federal  Highway  Administration  (Robson  1990) . 

4.  Care  must  be  taken  in  selecting  appropriate  positions  for  strain 
comparisons.  The  locations  for  which  strain  is  to  be  compared  are  essentially 
dependent  on  the  information  sought,  or  the  variables  that  are  considered  to 
be  unknown.  In  general,  the  number  of  strain  values  should  be  equal  to  or 
greater  than  the  number  of  unknown  quantities.  Locations  should  be  selected 
so  that  strain  responses  at  that  point  will  be  directly  affected  by  the  value 
of  an  unknown  parameter.  Assorted  load  configurations  can  also  be  used  to 
obtain  additional  strain  comparisons. 

5.  When  optimizing  various  structural  parameters,  it  is  important  to 
avoid  optimizing  two  or  more  variables  that  are  directly  correlated  or  have 
identical  effects  on  the  structural  response.  For  example,  the  flexural 
stiffness  of  a  girder  is  dependent  on  the  product  of  the  material  modulus  of 
elasticity  E  and  the  moment  of  inertia  I.  These  two  parameters  cannot  be 
optimized  at  the  same  time  if  stiffness  is  of  concern  since  they  have 
identical  effects  on  the  response.  While  attempting  to  optimize  the  flexural 
stiffness  of  a  beam,  one  term  might  approach  infinity  and  the  other  term  might 
approach  zero.  Another  important  factor  in  the  optimization  process  is  the 
constraints  on  a  given  parameter.  The  upper  and  lower  limits  should  be 
defined  using  good  engineering  judgment.  If  a  parameter  is  to  be  optimized  it 
is  necessary  that  the  range  be  sufficiently  broad  so  that  it  will  have  a 
measurable  effect  on  the  response.  However,  since  the  step  size  of  a 
parameter  change  is  based  on  a  percentage  of  the  total  range,  it  is  necessary 
that  the  variable  have  realistic  boundaries. 

6.  The  ability  of  the  program  to  obtain  a  good  correlation  is  dependent 
on  the  quality  of  the  model  and  the  relevance  of  the  unknown  parameters  to  the 
measured  response.  The  number  of  iterations  required  to  reach  the  desired 
correlation  is  directly  proportional  to  the  number  of  parameters  being  opti¬ 
mized.  It  is  desirable  to  keep  the  selected  unknown  variables  to  a  minimum  as 
it  reduces  the  computer  run  time  and  it  generally  improves  the  correlation. 
Fewer  strain  comparison  locations  are  required  if  the  number  of  optimized 
variables  is  reduced. 

7.  This  method  of  incorporating  an  optimization  algorithm  within  a 
structural  analysis  program  to  evaluate  structural  models  is  a  relatively  new 
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concept.  In  the  current  version  of  SAC,  a  rather  simplistic  process  was 
implemented  to  determine  if  the  concept  was  feasible.  Thus  far,  the  approach 
has  achieved  good  results,  so  it  may  be  desirable  to  adopt  a  more  sophisti¬ 
cated  optimization  algorithm.  More  complex  minimization  techniques  exist  that 
may  improve  the  convergence  time  and  correlation  of  the  models  with  their 
respective  experimental  data. 
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APPENDIX  B:  COMPUTER  IMPLEMENTATION 


1.  Even  though  it  is  desirable  to  keep  miter  gate  models  simple,  it  is 
possible  that  some  numeric  representations  will  be  relatively  large  compared  to 
most  PC  applications.  To  maintain  the  use  of  SAC  on  a  PC,  several  features  have 
been  incorporated  to  minimize  memory  requirements  and  maximize  performance  in 
speed  and  problem  size. 

2.  An  efficient  means  of  storing  data  arrays  has  been  implemented  in  SAC. 
A  dynamic  data  storage  scheme  eliminates  the  use  of  oversized  arrays  in  the 
program  routines  so  computer  memory  is  not  used  to  store  unnecessary  rows  and 
columns  of  zeros. 

3.  Typically,  with  large  finite  element  problems,  a  substantial  amount  of 
computer  memory  is  required  to  store  the  global  stiffness  matrix  which  is 
generated  in  the  execution  of  the  analysis  program.  This  problem  has  been 
minimized  in  SAC  with  the  utilization  of  an  efficient  skyline  solver.  The  intent 
of  a  skyline  solver  is  to  take  advantage  of  the  sparse  population  of  relevant 
terms  and  the  symmetry  of  the  stiffness  matrix.  Only  the  nonzero  terms  in  the 
upper  half  of  the  stiffness  matrix  are  stored.  Another  feature  employed  to 
reduce  memory  requirements  is  a  column  height  minimization  process.  This 
procedure  enhances  the  efficiency  of  the  skyline  solver  by  optimizing  the 
equation  numbers  associated  with  each  degree  of  freedom.  The  effect  is  that 
nonzero  column  heights  in  the  stiffness  matrix  are  minimized.  Therefore,  no 
effort  is  required  in  producing  meshes  with  efficient  nodal  numbering. 

4.  Perhaps  the  most  significant  feature  of  SAC  implementation  is  the  use 
of  the  Lahey  EM/32  compiler  and  linker.  This  feature  enables  the  program  to  use 
extended  memory  as  well  as  virtual  memory  for  data  storage,  thus  memory 
limitations  are  that  of  the  machine  rather  than  the  normal  amount  accessible  by 
a  DOS  operating  system.  This  compiler  also  takes  advantage  of  32-bit  processing 
of  386  and  486  computers  so  that  processing  speed  is  approximately  twice  that  of 
programs  created  with  standard  PC  compilers.  Utilization  of  the  Lahey  compiler 
and  operating  system  has  increased  the  level  of  PC  performance  tc  that  of  many 
workstations.  The  use  of  this  compiler,  however,  requires  that  the  program  be 
run  on  a  386  or  486  computer  with  sufficient  extended  memory  of  at  least 
1  megabyte.  SAC  is  extremely  portable  since  all  of  the  source  code  is  written 
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in  standard  FORTRAN  77.  The  program  can  be  loaded  on  any  main  frame  or  work¬ 
station  having  a  FORTRAN  77  compiler  with  minimal  effort. 
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APPENDIX  C:  ANALYSIS  INPUT  AND  RESULT  DATA 


1.  A  strong  effort  was  made  to  make  SAC  a  user-friendly  finite  element 
program.  Data  are  entered  into  SAC  with  the  use  of  an  input  file  of  ASCII 
format.  Input  data  are  grouped  into  individual  blocks  to  simplify  data  entry  and 
to  increase  organization.  The  data  blocks  can  be  presented  in  any  order  and  the 
data  in  each  block  are  entered  in  free  format.  Thus,  there  is  no  need  to  count 
columns  when  creating  or  editing  input  data.  The  entire  input  file  and  a 
condensed  result  file  for  the  John  Hollis  Bankhead  miter  gate  example  are  listed 
below.  These  file  listings  illustrate  the  input  requirements  and  flexibility  of 
SAC  as  well  as  the  output  presentation.  A  detailed  description  of  user  input  is 
presented  in  the  SAC  users  manual. 
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188 

329 

0 

38 

35 

36 

45 

330 

0 

38 

36 

37 

46 

331 

0 

37 

38 

39 

48 

332 

0 

37 

39 

40 

49 

333 

1 

30 

184 

185 

190 

336 

0 

30 

187 

188 

193 

337 

0 

37 

44 

45 

54 

338 

0 

37 

45 

46 

55 

339 

0 

36 

47 

48 

57 

340 

0 

36 

48 

49 

58 

341 

1 

30 

189 

190 

195 

344 

0 

30 

192 

193 

198 

345 

0 

36 

53 

54 

63 

346 

0 

36 

54 

55 

64 

347 

0 

36 

56 

57 

66 

348 

0 

36 

57 

58 

67 

349 

1 

30 

194 

195 

200 

352 

0 

30 

197 

198 

203 

353 

0 

36 

62 

63 

72 

354 

0 

36 

63 

64 

73 

Quoin  G.  C. 

Quoin  G.  B. 

Quoin  G.  A. 

Diaphragm  1  B 
Diaphragm  1  A 

Diaphragm  2  B 
Diaphragm  2  A 

Diaphragm  3  B 
Diaphragm  3  A 

Diaphragm  4  B 
Diaphragm  4  A 

Diaphragm  5  B 
Diaphragm  5  A 

Miter  G.  C. 

Miter  G.  B. 

Miter  G.  A. 

Miter  seal 

11  Quoin  V.  Plate  17-18 

12 

169  Skin  Plate  17-18 

172 

17  Miter  V.  Plate  17-18 

18 

20  Quoin  V.  Plate  16-17 

21 

174  Skin  Plate  16-17 

177 

26  Miter  V.  Plate  16-17 

27 

29  Quoin  V.  Plate  15-16 

30 

179  Skin  Plate  15-16 

182 

35  Miter  V.  Plate  15-16 

36 

38  Quoin  V.  Plate  14-15 

39 

184  Skin  Plate  14-15 

187 

44  Miter  V.  Plate  14-15 

45 

47  Quoin  V,  Plate  13-14 

48 

189  Skin  Plate  13-14 

192 

53  Miter  V.  Plate  13-14 

54 

56  Quoin  V.  Plate  12-13 

57 

194  Skin  Plate  12-13 

197 

62  Miter  V.  Plate  12-13 

63 

65  Quoin  V.  Plate  11-12 

66 

199  Skin  Plate  11-12 

202 

71  Miter  V.  Plate  11-12 

72 
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355 

0 

36 

65 

66 

75 

356 

0 

36 

66 

67 

76 

357 

1 

30 

199 

200 

205 

360 

0 

30 

202 

203 

208 

361 

0 

36 

71 

72 

81 

362 

0 

36 

72 

73 

82 

363 

0 

36 

74 

75 

84 

364 

0 

36 

75 

76 

85 

365 

1 

30 

204 

205 

210 

:>68 

0 

30 

207 

208 

213 

369 

0 

36 

80 

81 

90 

370 

0 

36 

81 

82 

91 

371 

0 

35 

83 

84 

93 

372 

0 

35 

84 

85 

94 

373 

1 

29 

209 

210 

215 

376 

0 

29 

212 

213 

218 

377 

0 

35 

89 

90 

99 

378 

0 

35 

90 

91 

100 

379 

0 

35 

92 

93 

102 

380 

0 

35 

93 

94 

103 

381 

1 

28 

214 

215 

220 

384 

0 

28 

217 

218 

223 

385 

0 

35 

98 

99 

108 

386 

0 

35 

99 

100 

109 

387 

0 

35 

101 

102 

111 

388 

0 

35 

102 

103 

112 

389 

1 

28 

219 

220 

225 

392 

0 

28 

222 

223 

228 

393 

0 

35 

107 

108 

117 

394 

0 

35 

108 

109 

118 

395 

0 

34 

110 

111 

120 

396 

0 

34 

111 

112 

121 

397 

1 

27 

224 

225 

230 

400 

0 

27 

227 

228 

233 

401 

0 

34 

116 

117 

126 

402 

0 

34 

117 

118 

127 

403 

0 

33 

119 

120 

129 

404 

0 

33 

120 

121 

130 

405 

1 

26 

229 

230 

235 

408 

0 

26 

232 

233 

238 

409 

0 

33 

125 

126 

135 

410 

0 

33 

126 

127 

136 

411 

0 

32 

128 

129 

138 

412 

0 

32 

129 

130 

139 

413 

1 

26 

234 

235 

240 

416 

0 

26 

237 

238 

243 

417 

0 

32 

134 

135 

144 

418 

0 

32 

135 

136 

145 

419 

0 

32 

137 

138 

147 

420 

0 

32 

138 

139 

148 

421 

1 

26 

239 

240 

245 

424 

0 

26 

242 

243 

248 

425 

0 

32 

143 

144 

153 

426 

0 

32 

144 

145 

154 

427 

0 

32 

146 

147 

156 

428 

0 

32 

147 

148 

157 

429 

1 

26 

244 

245 

250 

432 

0 

26 

247 

248 

253 

433 

0 

32 

152 

153 

162 

434 

0 

32 

153 

154 

163 

435 

1 

39 

4 

164 

439 

0 

39 

8 

168 

440 

1 

39 

13 

169 

444 

0 

39 

17 

173 

445 

1 

39 

22 

174 

449 

0 

39 

26 

178 

450 

1 

39 

31 

179 

454 

0 

39 

35 

183 

455 

1 

39 

40 

184 

74  Quoin  V.  Plate  10-11 

75 

204  Skin  Plate  10-11 

207 

80  Miter  V.  Plate  10-11 

81 

83  Quoin  V.  Plate  9-10 

84 

209  Skin  Plate  9-10 

212 

89  Miter  V.  Plate  9-10 

90 

92  Quoin  V.  Plate  8-9 

93 

214  Skin  Plate  8-9 

217 

98  Miter  V.  Plate  ^-V 

99 

101  Quoin  V.  Plate  7-8 

102 

219  Skin  Plate  7-8 

222 

107  Miter  V.  :  late  7-£ 

108 

110  Quoin  V.  Plate  6-7 

111 

224  Skin  Plate  6-7 

227 

116  Miter  V.  Plate  6-7 

117 

119  Quoin  V.  Plate  5-6 

120 

229  Skin  Plate  5-6 

232 

125  Miter  V.  Plate  5-6 

126 

128  Quoin  V.  Plate  4-5 

129 

234  Skin  Plate  4-5 

237 

134  Miter  V.  Plate  4-5 

135 

137  Quoin  V.  Plate  3-4 

138 

239  Skin  Plate  3-4 

242 

143  Miter  V.  Plate  3-4 

144 

146  Quoin  V.  Plate  2-3 

147 

244  Skin  Plate  2-3 

247 

152  Miter  V.  Plate  2-3 

153 

155  Quoin  V.  Plate  1-2 

156 

249  Skin  Plate  1-2 

252 

161  Miter  V.  Plate  1-2 

162 

Rigid  Links  to  Plates 
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459  0  39  44  188 

460  1  39  49  189 

464  0  39  53  193 

465  1  39  58  194 

469  0  39  62  198 

470  1  39  67  199 

474  0  39  71  203 

475  1  39  76  204 

479  0  39  80  208 

480  1  39  85  209 

484  0  39  89  213 

485  1  39  94  214 

489  0  39  98  218 

490  1  39  103  219 

494  0  39  107  223 

495  1  39  112  224 

499  0  39  116  228 

500  1  39  121  229 

504  0  39  125  233 

505  1  39  130  234 

509  0  39  134  238 

510  1  39  139  239 

514  0  39  143  243 

515  1  39  148  244 

519  0  39  152  248 

520  1  39  157  249 

524  0  39  161  253 

525  0  40  10  254  MITER  INCLINED  SUPPORT  ELEMENTS 

526  0  40  19  255 

527  0  40  28  256 

528  0  40  37  257 

529  0  40  46  258 

530  0  40  55  259 

531  0  40  64  260 

532  0  40  73  261 

533  0  40  82  262 

534  0  40  91  263 

535  0  40  100  264 

536  0  40  109  265 

537  0  40  118  266 

538  0  40  127  267 

539  0  40  136  268 

540  0  40  145  269 

541  0  40  154  270 

542  0  40  163  271 

543  0  41  272  275  DIAGONAL  N4-N161 

544  0  42  273  274  DIAGONAL  N8-N157 

545  0  39  4  272  RIGID  LINK  TO  DIAGONAL 

546  0  39  8  273  RIGID  LINK  TO  DIAGONAL 

547  0  39  157  274  RIGID  LINK  TO  DIAGONAL 

548  0  39  161  275  RIGID  LINK  TO  DIAGONAL 

OUTB  (SETUP  1:  CHAN  1-10,13-32) 


54 

42.0 

26.5 

G12  :  USF  -  MS 

54 

42.0 

-46.2 

G12  :  DSF  -  MS 

70 

42.0 

26.5 

G10  :  USF  -  MS 

70 

42.0 

-46.2 

G10  :  DSF  -  MS 

86 

42.0 

27.1 

G  8  :  USF  -  MS 

86 

42.0 

-45.4 

G  8  :  DSF  -  MS 

94 

42.0 

27.1 

G  7  :  USF  -  MS 

94 

42.0 

-45.4 

G  7  :  DSF  -  MS 

94 

42.0 

-8.4 

G  7  :  STIFFENER 

94 

42.0 

-45.4 

G  7  :  DSF  -  MS 

544 

580.0 

0.0 

X-BRACE  2 

544 

580.0 

0.0 

X- BRACE  2 

543 

580.0 

0.0 

X-BRACE  1 

543 

580.0 

0.0 

X-BRACE  1 

110 

42.0 

30.5 

G  5  :  USF  -  MS 

110 

42.0 

-42.0 

G  5  :  DSF  -  MS 

CIO 


226 

45.0 

33.2 

D5-6: 

USF 

- 

MS 

226 

45.0 

-39.3 

D5-6: 

DSF 

- 

MS 

126 

42.0 

32.1 

G  3  : 

USF 

- 

MS 

126 

42.0 

-40.4 

G  3  : 

DSF 

- 

MS 

142 

42.0 

33.6 

G  1  : 

USF 

- 

MS 

142 

42.0 

-40.0 

G  1  : 

DSF 

- 

MS 

87 

99.8 

38.5 

G  8  : 

USF 

- 

MITER  END 

87 

99.8 

-34.8 

G  8  : 

DSF 

- 

MITER  END 

95 

99.8 

38.5 

6  7  ; 

USF 

- 

MITER  END 

95 

99.8 

-34.8 

G  7  : 

DSF 

- 

MITER  END 

84 

42.0 

38.5 

G  8  : 

USF 

- 

QUOIN  END 

84 

42.0 

-34.8 

G  8  : 

DSF 

- 

QUOIN  END 

92 

42.0 

38.5 

G  7  : 

USF 

- 

QUOIN  END 

92 

42.0 

-34.8 

G  7  : 

DSF 

- 

QUOIN  END 

LOAD 

CASE  1  POOL  LEVEL  a  LIFT  LINE  2  -  CLICK  1 

WATER  LOAD 

0.0  229.6  157.6  3.611e-05  2 

ENDCASE  1 

CASE  2  POOL  LEVEL  a  LIFT  LINE  4  -  CLICK  3 

WATER  LOAD 

0.0  349.2  157.6  3.611e-05  2 

ENDCASE  2 

CASE  3  POOL  LEVEL  a  LIFT  LINE  6  -  CLICK  5 

WATER  LOAD 

0.0  469.2  157.6  3.611e-05  2 

ENDCASE  3 

CASE  4  POOL  LEVEL  a  LIFT  LINE  8  -  CLICK  7 

WATER  LOAD 

0.0  589,2  157,6  3.611e-05  2 

ENDCASE  4 

CASE  5  POOL  LEVEL  a  LIFT  LINE  10  -  CLICK  9 

WATER  LOAD 

0.0  709.2  157.6  3.611e-05  2 

ENDCASE  5 

CASE  6  POOL  LEVEL  3  LIFT  LINE  12  -  CLICK  11 

WATER  LOAD 

0.0  829.2  157.6  3.611e-05  2 

ENDCASE  6 

CASE  7  POOL  LEVEL  a  LIFT  LINE  14  -  CLICK  13 

WATER  LOAD 

0.0  949,2  157.6  3.611e-05  2 

ENDCASE  7 

#CASE  1  DEAD  LOAD  -  UNMITERED  POSITION 
GRAVITY  LOAD 
-1.0  2 
ENDCASE  1 
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Listing  of  condensed  result  file: 


MICRO  FINITE  ELEMENT  PROGRAM  -  SAC  PAGE  1 

John  Hollis  Bankhead  lock  and  dam  -  lower  miter  leaf  -  Black  Warrior  River,  Alab 


* INPUT  CONTROL  PARAMETERS: 


NUMBER  OF  NOOAL  POINTS  .  253 

NUMBER  OF  ELEMENTS . 524 

NUMBER  OF  MATERIAL  GROUPS  .  39 

SPATIAL  DIMENSIONS  .  3 

MAX  DEGREES  OF  FREEDOM/NODE  .  6 

MAX  NUMBER  OF  NOOES/ELEMENT  .  4 

NUMBER  OF  LOAD  CASES .  1 

NUMBER  OF  BEAMS  FOR  OUTPUT .  0 

NUMBER  OF  PLATES  FOR  OUTPUT  .  0 

NUMBER  OF  PARAMETERS  FOR  OPTIMIZATION  .  .  0 


•MATERIAL  PROPERTIES 


MATERIAL  GROUP  .  1 

2-NOOE  EULER-BERNOULLI  BEAM 

ELASTIC  MODULUS  (E)  .  2.900E>04 

SHEAR  MODULUS  (G)  .  1.200E+04 

MOMENT  OF  INERTIA  ABOUT  y  .  .  3.877E+04 
MOMENT  OF  INERTIA  ABOUT  z  .  .  5.000E'f03 
TORSIONAL  SHAPE  FACTOR  (J).  .  1.281E+01 
WARPING  TORSION  CONST.  (Cw)  .  1.853E4-05 

CROSS-SECTIONAL  AREA .  5.775E+01 

ROTATION  ABOUT  BEAM  AXES.  .  .  O.OOOE+00 
MATERIAL  UNIT  WEIGHT .  2.830E-04 

MATERIAL  GROUP  .  2 

2-NOOE  EULER-BERNOULLI  BEAM 

ELASTIC  MODULUS  (E)  .  2.900E+04 

SHEAR  MODULUS  (G)  .  1.200E>04 

MOMENT  OF  INERTIA  ABOUT  y  .  .  5.381E>04 
MOMENT  OF  INERTIA  ABOUT  z  .  .  5.000E+03 
TORSIONAL  SHAPE  FACTOR  (J).  .  2.026E+01 
WARPING  TORSION  CONST.  <Cw)  .  3.879E+05 

CROSS-SECTIONAL  AREA .  6.825E+01 

ROTATION  ABOUT  BEAM  AXES.  .  .  O.OOOE+00 
MATERIAL  UNIT  WEIGHT .  2.830E-04 


MATERIAL  GROUP  .  26 

4 -NODE  RECTANGULAR  KIRCHOFF  PLATE 

ELASTIC  MODULUS  .  1.000E-03 

POISSON  RATIO  .  1.500E-01 

ELEMENi  THICKNESS  .  3.750E-01 

MATERIAL  UNIT  WEIGHT .  2.830E-04 
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MICRO  FINITE  ELEMENT  PROGRAM  -  SAC  PAGE  11 

John  Hollis  Bankhead  lock  and  dam  -  lower  miter  leaf  -  Black  Warrior  River,  Alab 


*N00AL  COORDINATES 


NODE 

X- COORD 

Y-COORD 

2 -COORD 

1 

O-OODOOE-t-OO 

O.OOOOOE^OO 

0.00000E«00 

2 

2.17000E+01 

O.OOOOOE+00 

7.23300E+00 

3 

3.51000E+01 

O.OOOOOE+00 

1.17000E+01 

4 

8.65000E+01 

O.OOOOOE+00 

3.20000E«01 

5 

2.28250E+02 

0.00000E«00 

3.20000E+01 

6 

3.70000E+02 

O.OOOOOE-i^OO 

3.20000E+01 

7 

5.11750E+02 

O.OOOOOE+00 

3.20000E'i-01 

8 

6.53500E+02 

O.OOOOOE^^OO 

3.20000E'<^01 

9 

7.04900E+02 

O.OOOOOE^OO 

1.17000E+01 

10 

7.40000E+02 

O.OOOOOE+00 

O.OOOOOE+00 

11 

O.OOOOOE+00 

4.80000E>01 

O.OOOOOE+OO 

12 

3.51000E+01 

4.80000E+01 

1.17000E+01 

13 

8.65000E+01 

4.80000E+01 

3.20000E«01 

14 

2.28250E402 

4.80000E+01 

3.20000E+01 

15 

3.70000E+02 

4.80000E+01 

3.20000E>01 

16 

5.11750E+02 

4.80000E+01 

3.20000E>01 

17 

6.53500E+02 

4.80000E>01 

3.20000E«01 

18 

7.04900E+02 

4.80000E+01 

1.17000E+01 

19 

7.40000E+02 

4.80000E'»01 

O.OOOOOE+00 

20 

O.OOOOOE-^OO 

9.60000E+01 

O.OOOOOE^OO 

21 

3.51000E+01 

9.60000E+01 

1.17000E+01 

22 

8.65000E+01 

9.60000E+01 

3.20000E+01 

23 

2.28250E+02 

9.60000E'<’01 

3.20000E+01 

24 

3.70000E+02 

9.60000E+01 

3.20000E+01 

25 

5.11750E+02 

9.60000e-f01 

3.2aOQO£*Ol 

26 

6.53500E+02 

9.60000E+01 

3.20000E+01 

27 

7.04900E+02 

9.60000E-i’01 

1.17000E+01 

28 

7.40000E+02 

9.60000E+01 

O.OOOOOE-^OO 

29 

O.OOOOOE+00 

1.44000E+02 

0.00000E-»00 

30 

3.51000E+01 

1.44000E+02 

1.17000E+01 

31 

8.65000E+01 

1.44000E+02 

3.20000E+01 

32 

2.28250E+02 

1.44000E+02 

3.20000E+01 

33 

3.70000E+02 

1.44000E+02 

3.20000E+01 

34 

5.11750E+02 

1.44000E+02 

3.20000E+01 

35 

6.53500E+02 

1.44000E+02 

3.20000E+01 

36 

7.04900E+02 

1.44000E+02 

1.17000E+01 

37 

7.40000E+02 

1 .44000E+02 

O.OOOOOE+00 

38 

O.OOOOOE+00 

1 .92000E+02 

0.00000E>00 

39 

3.51000E+01 

1.92000E+02 

1.17000E+01 

40 

8.65000E+01 

1 .92000E+02 

3.20000E+01 

41 

2.28250E+02 

1.92000E+02 

3.20000E+01 

42 

3.70000E+02 

1.92000E+02 

3.20000E+01 

43 

5.11750E+02 

1.92000E+02 

3.20000E+01 

44 

6.53500E+02 

1.92000E+02 

3.20000E+01 

45 

7.04900E+02 

1.92000E+02 

1.17000E+01 

46 

7.40000E+02 

1 .92000E+02 

O.OOOOOE^OO 

47 

O.OOOOOE-t-00 

2.40000E+02 

O.OOOOOE^OO 
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♦ELEMENT  TOPOLOGY 


EMENT 

MATERIAL 

1  NO 

2  NO 

3  NO 

4  NO 

1 

13 

1 

2 

2 

13 

2 

3 

3 

13 

3 

4 

4 

14 

4 

5 

5 

15 

5 

6 

6 

15 

6 

7 

7 

14 

7 

8 

8 

13 

8 

9 

9 

13 

9 

10 

10 

10 

11 

12 

11 

10 

12 

13 

12 

11 

13 

14 

13 

12 

14 

15 

14 

12 

15 

16 

15 

11 

16 

17 

16 

10 

17 

18 

17 

10 

18 

19 

18 

10 

20 

21 

19 

10 

21 

22 

20 

11 

22 

23 

299 

38 

1 

3 

12 

11 

300 

38 

3 

4 

13 

12 

301 

31 

164 

165 

170 

169 

302 

31 

165 

166 

171 

170 

303 

31 

166 

167 

172 

171 

304 

31 

167 

168 

173 

172 

305 

38 

8 

9 

18 

17 

306 

38 

9 

10 

19 

18 

307 

38 

11 

12 

21 

20 

308 

38 

12 

13 

22 

21 

309 

31 

169 

170 

175 

174 

310 

31 

170 

171 

176 

175 

311 

31 

171 

172 

177 

176 

312 

31 

172 

173 

178 

177 

313 

38 

17 

18 

27 

26 

314 

38 

18 

19 

28 

27 

315 

38 

20 

21 

30 

29 

316 

38 

21 

22 

31 

30 

317 

31 

174 

175 

180 

179 

318 

31 

175 

176 

181 

180 

319 

31 

176 

177 

182 

181 

320 

31 

177 

178 

183 

182 

321 

38 

26 

27 

36 

35 

322 

38 

27 

28 

37 

36 

323 

38 

29 

30 

39 

38 

324 

38 

30 

31 

40 

39 

325 

31 

179 

180 

185 

184 

326 

31 

180 

181 

186 

185 

327 

31 

181 

182 

187 

186 

328 

31 

182 

183 

188 

187 

329 

38 

35 

36 

45 

44 

330 

38 

36 

37 

46 
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*NODAL  RESTRAINT  CONDITIONS 


NODE 


X-DISP  Y-DISP  Z-DISP  X-ROT 


Y-ROT  Z-ROT 


2  1 

7  0 

155  1 


1  1 
0  1 
0  1 


0  0  0 
0  0  0 
1  0  1 


‘SKYLINE  OPTIMIZATION  RESULTS 
INITIAL  VALUES; 

MAXIMUM  COLUMN  HEIGHT:  961 

SIZE  OF  GLOBAL  STIFFNESS  VECTOR:  3708.432  KB 

OPTIMIZED  VALUES: 

MAXIMUM  COLUMN  HEIGHT:  143 

SIZE  OF  GLOBAL  STIFFNESS  VECTOR:  1098.904  KB 
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*******  load  data  for  load  case  #  1 


‘APPLIED  GRAVITY  LOADS 


GRAVITY  FACTOR  .  -1.00 

GRAVITY  DIRECTION  CODE  .  2 


***********  OUTPUT  FOR  LOAD  CASE  #  1  *********** 
•••••*******«•****•**••••••••******•*•***••••••* 


‘NODAL  DISPLACEMENTS  IN  GLOBAL  COORDINATES 


NODE 

X-DlSP 

Y-DISP 

Z-DISP 

X-ROT 

Y-ROT 

Z-ROT 

1 

0.1050E-02 

-.3308E-02 

-.2263E-02 

-.2430E-03 

-.1068E-03 

0.1298E-03 

2 

O.OOOOE+00 

0.000OE‘OO 

0.0000E‘00 

-.1843E-03 

-.1112E-03 

-.1278E-03 

3 

-.1774E-02 

-.3069E-02 

0.1096E-02 

-.1139E-03 

-.1151E-03 

-.3582E-03 

4 

-.7381E-02 

-.2925E-01 

0.5657E-02 

-.5105E-03 

-.1058E-03 

-.6095E-03 

5 

-.1202E-01 

-.1333E+00 

0.1405E-01 

-.1447E-02 

-.1427E-04 

-.5573E-03 

6 

-.1317E-01 

-.2178E‘00 

0.1128E-01 

-.2273E-02 

0.5214E-04 

-.4978E-03 

7 

-.1300E-01 

-.2888E+00 

0.0000E‘00 

-.3047E-02 

0.1069E-03 

-.4445E-03 

8 

-.1242E-01 

-.3486E+00 

-.2030E-01 

-.3763E-02 

0.1897E-03 

-.4013E-03 

9 

-.1640E-01 

-.4453E+00 

-.3077E-01 

-.3983E-02 

0.2133E-03 

-.3055E-03 

10 

-.1893E-01 

-.5022E+00 

-.3839E-01 

-.4063E-02 

0.2190E-03 

-.2584E-03 

11 

0.4811E-02 

-.9313E-02 

0.5650E-03 

-.5145E-04 

0.2312E-03 

-.1316E-03 

12 

0.7995E-02 

-.1395E-01 

-.7392E-02 

-.1712E-03 

0.2312E-03 

-.2232E-03 

13 

0.1251E-01 

-.3020E-01 

-.1942E-01 

-.5272E-03 

0.2350E-03 

-.4345E-03 

14 

0.9933E-02 

-.1334E+00 

-.5542E-01 

-.1447E-02 

0.2767E-03 

-.4794E-03 

15 

0.8387E-02 

-.2178E+00 

-.9784E-01 

-.2274E-02 

0.3211E-03 

-.4524E-03 

16 

0.7540E-02 

-.2887E+00 

-.1463E+00 

-.3049E-02 

0.3621E-03 

-.4257E-03 

17 

0.7086E-02 

-.3486E+00 

-.2009E‘00 

-.3767E-02 

0.4090E-03 

-.3969E-03 

18 

-.1376E-02 

-.4453E+00 

-.2222E+00 

-.3999E-02 

0.4172E-03 

-.3037E-03 

19 

-.6292E-02 

-.5021E+00 

-.2368E‘00 

-.4073E-02 

0.4187E-03 

-.2467E-03 

20 

0.1174E-01 

-.1551E-01 

0.1285E-02 

-.1204E-03 

0.5454E-03 

-.1520E-03 

21 

0.1814E-01 

-.2087E-01 

-.1786E-01 

-.2279E-03 

0.5458E-03 

-.2512E-03 

22 

0.2925E-01 

-.3329E-01 

-.4596E-01 

-.5676E-03 

0.5482E-03 

-.4087E-03 

23 

0.2832E-01 

-.1337E+00 

-.1249E+00 

-.1450E-02 

0.5677E-03 

-.4410E-03 
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*  INTERNAL  FORCES  FOR  BEAMS 


ELEMENT 

NODE 

AXIAL 

SHEAR -Y 

SHEAR -Z 

TORSION 

HOHENT-Y 

MOMENT -Z 

1 

1 

27.70 

-253.05 

26.36 

437.83 

50.07 

-1226.97 

2 

-27.70 

253.05 

-26.36 

-437.83 

-653.06 

-4561.11 

2 

2 

213.98 

295.73 

-18.26 

438.00 

653.06 

4561.10 

3 

-213.98 

-295.73 

18.26 

-438.00 

-395.11 

-383.94 

3 

3 

144.88 

78.87 

-26.19 

557.89 

408.61 

2410.28 

4 

-144.88 

-78.87 

26.19 

-557.89 

1038.86 

1948.19 

4 

4 

63.12 

13.05 

0.66 

43.49 

-935.71 

871.70 

5 

-63.12 

-13.05 

-0.66 

-43.49 

841.88 

978.66 

5 

5 

17.07 

5.96 

0.53 

61.53 

-753.79 

361.70 

6 

-17.07 

-5.96 

-0.53 

-61.53 

679.00 

483.43 

6 

6 

-2.51 

2.52 

0.10 

57.65 

-597.58 

123.87 

7 

2.51 

-2.52 

-0.10 

-57.65 

583.40 

232.86 

7 

7 

-7.90 

-0.04 

-4.22 

33.23 

-506.13 

-46.68 

8 

7.90 

0.04 

4.22 

-33.23 

1103.65 

41.69 

8 

8 

-6.02 

-1.34 

8.78 

291.04 

-1037.27 

-57.91 

9 

6.02 

1.34 

-8.78 

-291.04 

551.94 

-15.97 

9 

9 

-0.39 

-2.92 

13.90 

362.92 

-543.99 

-130.62 

10 

0.39 

2.92 

-13.90 

-362.92 

29.57 

22.76 

10 

11 

-35.61 

-9.72 

0.24 

1359.01 

-3.54 

12.28 

12 

35.61 

9,72 

-0.24 

-1359.01 

-5.19 

-371.86 

11 

12 

10.25 

66.28 

-5.81 

1170.72 

3.95 

2003.93 

13 

-10.25 

-66.28 

5.81 

-1170.72 

317.17 

1658.89 

12 

13 

37.85 

23.45 

-1.74 

66.81 

-326.07 

1708.07 

14 

-37.85 

-23.45 

1.74 

-66.81 

572.09 

1616.15 

13 

14 

25.55 

11.26 

0.40 

106.95 

-572.65 

770.54 

15 

-25.55 

-11.26 

-0.40 

-106.95 

516.34 

825.89 

14 

15 

13.99 

5.29 

0.18 

100.27 

-516,31 

347.49 

16 

-13.*^ 

-5.29 

-0.18 

-100.27 

490.50 

402.06 

15 

16 

6.68 

0.97 

-0.18 

52.18 

-490.94 

38.96 

17 

-6.68 

-0.97 

0.18 

-52.18 

516.63 

97.92 

16 

17 

3.14 

-1.61 

6.11 

716.75 

-513.66 

-48.16 

18 

-3.14 

1.61 

-6.11 

-716.75 

176.19 

-40.75 

17 

18 

1.72 

-3.04 

4.64 

838.37 

-175.36 

-176.82 

19 

-1.72 

3.04 

-4.64 

-838.37 

3.65 

64.17 

18 

20 

-1.52 

11.14 

-1.22 

1269.92 

-2.86 

441.55 

21 

1.52 

-11.14 

1.22 

-1269.92 

48.13 

-29.40 

19 

21 

-0.34 

36.46 

-1.86 

1071.04 

-48.90 

1064.53 

22 

0.34 

-36.46 

1.86 

-1071.04 

151.71 

950.61 

20 

22 

13.69 

24.52 

-0.73 

64.08 

-157.27 

1771.13 

23 

-13.69 

-24.52 

0.73 

-64.08 

260.74 

1704.92 

21 

23 

15.37 

14.06 

-0.49 

106.90 

-262.18 

976.85 

24 

-15.37 

-14,06 

0.49 

-106.90 

331.20 

1016.46 

22 

24 

11.93 

7.48 

0.22 

100.83 

-330.74 

513.03 

25 

-11.93 

-7.48 

-0.22 

-100.83 

299.30 

547.79 

23 

25 

7.70 

2.31 

0.88 

53.07 

-298.63 

144.92 

26 

-7.70 

-2.31 

-0.88 

-53.07 

174.42 

182.43 

24 

26 

3.92 

0.28 

1.86 

TOO. 20 

-170.77 

6.52 

27 

-3.92 

-0.28 

-1.86 

-780.20 

67.71 

8.84 

25 

27 

1.60 

-1.29 

1.75 

934.70 

-67.11 

-152.63 

28 

-1.60 

1.29 

-1.75 

-934.70 

2.36 

105.00 

26 

29 

0.13 

15.38 

-1.31 

1218.83 

-2.06 

482.34 

30 

-0.13 

-15.38 

1.31 

-1218.83 

50.69 

86.55 

27 

30 

2.95 

25.72 

-1.54 

1011.51 

-51.22 

773.17 

31 

-2.95 

-25.72 

1.54 

-1011.51 

136.23 

647.96 

28 

31 

7.24 

23.53 

0.33 

61.82 

-139.97 

1693.63 

32 

-7.24 

-23.53 

-0.33 

-61.82 

93.18 

1641.61 
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•SUMMATION  OF 

APPLIED  FORCES 

FIX) 

F(Y) 

F(Z) 

M(X) 

M(Y) 

M(Z) 

O.OOOE+00 

-5.492E+02 

O.OOOE+00 

O.OOOE+OO 

0.000E*00 

O.OOOE^OO 

•NOOAL  REACTIONS  AT  FIXED 

NODES 

NODE  # 

F(X) 

F(Y) 

F(Z) 

M(X) 

M(Y)  M(Z) 

2 

190.831 

549.183 

16.579 

0.000 

0.000  0.000 

7 

0.000 

0.000 

-4.315 

0.000 

0.000  0.000 

155 

-190.831 

0.000 

-12.264 

-609.411 

0.000  -2203.841 

SUMMATION 

OF  APPLIED 

FORCES  AND  REACTIONS 

FIX) 

F(Y) 

F<Z) 

M<X) 

M(Y) 

M(Z) 

-4.369E- 

06  -5.511E- 

07  -2.448E- 

04  -6.094E- 

08  O.OOOE+00 

-2.204E-06 

TOTAL  RUN  TIME  (HH;MM:SS.hh)  =  00:02:55.98 
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APPENDIX  D:  A  DEVELOPMENT  OF  WARPING  TORSION  TERMS 


1.  Restraint  of  torsional  flange  warping  can  be  represented  in  two  ways. 
One  method  is  to  use  a  three-dimensional  modeling  approach  in  which  the  member 
flanges  are  modeled  independently.  This  approach  results  in  independent  dis¬ 
placements  of  each  flange  rather  than  having  only  a  single  set  of  displacements 
for  a  given  girder  or  diaphragm  cross  section.  The  other  approach  is  to 
implement  restrained  warping  torsion  terms  into  the  torsional  stiffness  of  the 
frame  elements.  The  development  of  warping  torsion  stiffness  terms  is  available 
in  most  advanced  steel  design  and  structural  analysis  references.  However, 
warping  constants  are  not  often  incorporated  in  ordinary  frame  analysis  since  the 
connection  details  of  most  structures  do  not  warrant  consideration  of  warping 
torsion,  and  its  effect  is  generally  not  significant  to  overall  frame  behavior. 
Following  is  a  brief  development  of  the  warping  torsion  terms  that  may  be  applied 
in  the  element  stiffness  matrices  for  the  examples  which  consider  warping 
torsion. 

2.  The  total  torsional  resistance  of  a  given  I-shaped  cross  section  is  the 
summation  of  St.  Venant's  torsion  and  the  restrained  flange  warping  torsion. 

r  =  Tj  +  r2  (Di) 

where  is  St.  Venant's  torsion  and  T2  is  the  restrained  warping  torsion. 

The  expression  for  St.  Venant's  torsion  is 


where  G  is  the  shear  modulus,  J  the  polar  moment  of  inertia,  0  the  angle  of 
twist,  and  L  the  torsional  unbraced  length. 

3.  The  warping  torsion  terms  are  due  to  the  bending  resistance  of  the 
flanges  about  their  strong  axes.  The  forces  that  are  generated  when  the  flanges 
are  braced  against  warping  are  illustrated  in  Figure  D1 .  Applying  basic  beam 
theory  to  obtain  the  resistance  of  a  single  flange  yields 


D1 


^1 

^1 

/ 

/ 

/ 

m 

^2! 

,  / 

‘  / 

/ 

,  y 

'  / 

^  C 

i 

L 

- ’^0 

_ r 

TORQUE  DUE  TO 

FLANGE  REACTIONS 

r  n 

REACTIONS  OF  FLANGE  1 

DUE  TO  RESTRAINED  WARPING 

Figure  Dl.  Restrained  warping  torsion  In  beam  flanges 


12El£  Aj 

P 


and 


6EI^  Af 

P 


(D3) 


where  the  subscript  f  designates  flange,  Is  the  flange  moment  of  Inertia  about 
Its  strong  axis,  Ff  and  Hf  represent  the  flange  force  and  moment  resultants  due 
to  a  displacement  Zi,  and  L  is  the  torsional  unbraced  length.  The  torque  applied 
by  the  warping  torsion  T2  is  then  the  summation  of  the  flange  forces  F^  times 
their  respective  distances  to  the  neutral  axis  . 

Substitution  of  the  expression  for  Fj  of  Equation  D3  into  Equation  D4  yields 


D2 


^2  =  Ei 


(D5) 


12  El,  A,  y, 

1? 


For  an  angle  of  twist  0  assuming  small  displacement  theory 


0  w  tan  0  =  _i 

/i 


(D6) 


Therefore , 


Ai  =  0 


(D7) 


Substituting  the  expression  for  into  Equation  D5  yields 


Tz 


12  E  e  X' 


71 


2 


(D8) 


The  summation  terns  in  Equation  D8  are  known  as  the  warping  constant  of  a  cross 
section 

=  7i" 


4.  For  symmetric  I  sections  it  is  assumed  that  the  moment  of  inertia  about 
the  weak  axis  ly  is  due  only  to  the  flanges,  simplifying  the  calculation  of  the 
warping  constant  as 

r  (DIO) 

4 


where  d  is  the  depth  of  the  I  section.  Substituting  Equations  D2 ,  D8  and  D9  into 
Equation  D1  yields  the  magnitude  of  the  overall  torsional  resistance  T. 
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User's  Guide:  Computer  Program  with  Interactive  Graphics  for 
Analysis  of  Plane  Frame  Structures  (CFRAME) 

User's  Guide;  Computer  Program  for  Generation  of  Engineering 
Geometry  (SKETCH) 

User's  Guide:  Computer  Program  to  Calculate  Shear,  Moment, 
and  Thrust  (CSMT)  from  Stress  Results  of  a  Two-Dimensional 
Finite  Element  Analysis 

Basic  Pile  Group  Behavior 

Reference  Manual;  Computer  Graphics  Program  for  Generation  of 
Engineering  Geometry  (SKETCH) 

Case  Study  of  Six  Major  General-Purpose  Finite  Element  Programs 

User’s  Guide:  Computer  Program  for  Optimum  Dynamic  Design 
of  Nonlinear  Metal  Plates  Under  Blast  Loading  (CSDOOR) 

User’s  Guide;  Computer  Program  for  Determining  Induced 
Stresses  and  Consolidation  Settlements  (CSETT) 

Seepage  Analysis  of  Confined  Row  Problems  by  the  Method  of 
Fragments  (CFRAG) 

User’s  Guide  for  Computer  Program  CGFAG,  Concrete  General 
Flexure  Analysis  with  Graphics 

Computer-Aided  Drafting  and  Design  for  Corps  Structural 
Engineers 

Decision  Logic  Table  Formulation  of  ACI 31 8-77,  Building  Code 
Requirements  for  Reinforced  Concrete  for  Automated  Con¬ 
straint  Processing,  Volumes  I  and  II 

A  Case  Committee  Study  of  Rnite  Element  Analysis  of  Concrete 
Flat  Slabs 

User’s  Guide:  Computer  Program  for  Two-Dimensional  Analysis 
of  U-Frame  Structures  (CUFRAM) 

User’s  Guide:  For  Concrete  Strength  Investigation  and  Design 
(CASTR)  in  Accordance  with  ACI  318-83 

Rnite-Element  Method  Package  for  Solving  Steady-State  Seepage 
Problems 

User's  Guide;  A  Three  Dimensional  Stability  Analysis/Design 
Program  (3DSAD)  Module 
Report  1 ;  Revision  1 :  General  Geometry 
Report  2:  General  Loads  Module 
Report  6:  Free-Body  Module 
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Iner':.  Guid?::  2  D  rramo  Analycis  Link  Pronram  (LINK2D) 


Finite  Liemcnt  Stiidies  ot  a  Horizontally  F-'ramed  Miter  Gate 

Report  1 ;  Initial  and  Refined  Finite  Element  Models  (Phases 
A,  B.  and  C),  Vclumos  I  and  il 
Report  ?;  Simplified  Frame  Model  (Phase  D) 

:!l  3;  Aitfirnato  Configuration  Miter  Gate  Finite  Element 
Studies-Open  Section 

Aitornato  Configuration  Miter  Gate  Finite  Element 
Studies-Closed  Sections 

Report  5:  Alternate  Configuration  Miter  Gate  Finite  Element 
Studies-Additional  Closed  Sections 
Report  6;  Elastic  Buckling  of  Girders  in  Horizontally  Framed 
Miter  Gates 

Report  7:  Application  and  Summary 


Repor 


Report  4; 


User's  Guide:  UTEXAS2  Slope-Stability  Package;  Volume  L 
User’s  Manual 


Instruclicn  Report  ITL-S7-5 
instruction  Report  ITL-87-6 


Technical  Report  !TL-87'8 


Sliding  Stability  of  Concrete  Structures  (CSLIDE) 

Criteria  Specifications  for  and  Validation  of  a  Computer  Program 
for  the  Design  or  Investigation  of  Horizontally  Framed  Miter 
Gates  (CMITER) 

Procedure  for  Static  Analysis  of  Gravihy  Dam.s  Using  the  Finite 
Element  Method  -  Phase  la 


Instruction  Report  iTL-88-1 
Tochnicol  Report  ITL-88-1 
Technical  Report  ITL-88-2 
Instruction  Report  iTL-88'2 


Instruction  Report  ITL-88-4 


Instruction  Report  GL-87-1 
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Technical  Report  ITL-8S 


User's  Guide:  Computer  Program  for  Analysis  of  Planar  Grid 
Structures  (CGRID) 

Development  of  Design  Formulas  for  Ribbed  Mat  Foundations 
on  Expansive  Soils 

User's  Guide:  Piie  Group  Graphics  Display  (CPGG)  Post¬ 
processor  to  CPGA  Program 

User's  Guide  for  Design  and  Investigation  of  Horizontally  Framed 
Miter  Gates  (CMITER) 

User’s  Guide  for  Revised  Computer  Program  to  Calculate  Shear, 
Moment,  and  Thrust  (CSMT) 

User's  Guide;  UTEXAS2  Slope-Stability  Package;  Volume  II, 
Theory 

User's  Guide:  Pile  Group  Analysis  (CPGA)  Computer  Group 

CBASIN-Structural  Design  of  Saint  Anthony  Falls  Stilling  Basins 
According  to  Corps  of  Engineers  Criteria  for  Hydraulic 
Structures:  Com.puter  Program  X0098 
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Instruction  Report  ITL-90-2 
Technical  Report  ITL-91-3 

Instruction  Report  ITL-91-1 

Instruction  Report  I7L-87-2 
(Revised) 
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Technical  Report  ITL-92-4 
Instruction  Report  ITL-92-3 

Instruction  Report  ITL-92-4 
Instruction  Report  ITL-92-5 
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CCHAN-Structural  Design  of  Rectangular  Channels  According 
to  Corps  of  Engineers  Criteria  for  Hydraulic 
Structures:  Computer  Program  X0097 

The  Response-Spectrum  Dynamic  Analysis  of  Gravity  Dams  Using 
the  Finite  Element  Method;  Phase  II 

State  of  the  Art  on  Expert  Systems  Applications  in  Design, 
Construction,  and  Maintenance  of  Structures 

User’s  Guide:  Computer  Program  for  Design  and  Analysis 
of  Sheet  Pile  Walls  by  Classical  Methods  (CWALSHT) 

Investigation  and  Design  of  U-Frame  Structures  Using 
Program  CUFRBC 

Volume  A:  Program  Criteria  and  Documentation 
Volume  B:  User’s  Guide  for  Basins 
Volume  C:  User’s  Guide  for  Channels 

User’s  Guide:  Computer  Program  for  Two-Dimensional  Analysis 
of  U-Frame  or  W-Frame  Structures  (CWFRAM) 

User’s  Guide:  Pile  Group-Concrete  Pile  Analysis  Program 
(CPGC)  Preprocessor  to  CPGA  Program 

Application  of  Rnite  Element,  Grid  Generation,  and  Scientific 
Visualization  Techniques  to  2-D  and  3-D  Seepage  and 
Groundwater  Modeling 

User’s  Guide:  Computer  Program  for  Design  and  Analysis 
of  Sheet-Pile  Walls  by  Classical  Methods  (CWALSHT) 

Including  Rowe’s  Moment  Reduction 

User’s  Guide  for  Concrete  Strength  Investigation  and  Design 
(CASTR)  in  Accordance  with  ACI  318-89 

Rinite  Element  Modeling  of  Welded  Thick  Plates  for  Bonneville 
Navigation  Lock 

Introduction  to  the  Computation  of  Response  Spectrum  for 
Earthquake  Loading 

Concept  Design  Example,  Computer  Aided  Structural 

Modeling  (CASM) 

Report  1 :  ^heme  A 
Report  2:  Scheme  B 
Report  3:  Scheme  C 

User’s  Guide:  Computer-Aided  Structural  Modeling 

(CASM)  -  Version  3.00 

Tutorial  Guide:  Computer-Aided  Structural  Modeling 

(CASM)  -  Version  3.00 
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Technical  Report  ITL-92-1 2 
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Optimization  of  Steel  Pile  Foundations  Using  Optimality  Criteria 

Refined  Stress  Analysis  of  Melvin  Price  Locks  and  Dam 

Knowledge-Based  Expert  System  for  Selection  and  Design 
of  Retaining  Structures 

Evaluation  of  Thermal  and  Incremental  Construction  Effects 
tor  Monoliths  AL-3  and  AL-5  of  the  Melvin  Price  Locks 
and  Dam 

User’s  Guide:  UTEXAS3  Slooe-Stability  Package;  Volume  IV, 

User's  Manual 

Computer-Aioed.  Field-Verified  Structural  Evaluation 

Report  1 :  Deveiooment  of  Comouter  Modeiinp  Techniques 
tor  Miter  Lock  Gates 
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Jun  1992 
Seo 1992 
Seo  1992 
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